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What is a Physicist? 


ECENTLY, at a meeting of responsible 

Army, Navy and business officers to con- 
sider problems of national defense, the need for 
physicists to help in these problems was dis- 
cussed. Several of those present spoke in no un- 
certain terms indicating that they had no need 
for men who dabbled only in atoms and who 
gave no thought to concrete problems whose 
immediate solution was required for the adequate 
defense of the country. This incident is men- 
tioned since it illustrates very well the reputation 
which physicists have acquired in the minds of 
many of the leaders of our country. 

The rather general belief that physicists are 
interested only in atoms probably has arisen 
because the progress in atomic and nuclear phys- 
ics is so spectacular that newspapers print news 
in this field to the almost complete exclusion of 
news of any other work of the physicist. As a 
result, the painstaking care which the physicist 
has used in the establishment of precise standards 
of temperature, in the analysis of internal fric- 
tion in vibrating rods, 
in the exact specifica- 


This situation is unfortunate and is distinctly 
harmful to the development of physics as a pro- 
fession. It often leads to the hiring of others to 
do the work that the physicist has been espe- 
cially trained to do. One doesn’t call in a dentist 
to treat pneumonia just because he has the title 
of doctor. Yet many an industrialist will call in 
any kind of technical man except the physicist 
to measure accurately the coefficient of expan- 
sion or to interpret an x-ray diffraction pattern 
or to carry out a spectroscopic analysis—jobs for 
which the physicist has more specific training 
than does the physician for the treatment of 
pneumonia. 

What is the physicist to do about this? 
First of all, he must realize that no one else 
will do anything unless he makes the first 
move. It has been suggested that the first move 
is to formulate a comprehensive definition of 
the work of the physicist and then to make a 
concerted effort to bring this definition to the 
attention of the right people. To define the 

work of the physicist 





tion of color, in the de- 
velopment of optical in- 
struments, in the design 
of electrical measuring 
instruments and in in- 
numerable other proj- 
ects of the greatest im- 
mediate importance has 
escaped the attention 
of those who are not 
directly concerned with 
the progress in physics. 


LEY S. BALLARD. 


spring. 





For December 


Qualitative Spectrochemical 
Agriculture and Geochemistry, by STAN- 


Also contributed research papers on: 


Flow of Air Through Porous Media, Ther- 
moplastic Behavior, Absorption Spec- 
tra, Theory of the Cyclotron, Noncubic 
Lattice Constants, Spark Gaps with 
Short Time Lag, and Isochronous Hair- 


is not easy. Your ed- 
itor will be pleased to 
receive suggestions for 
a concise but accu- 
rate definition of the 
physicist’s work. Per- 
haps there are better 
ways to emphasize 
the broad nature of 
a physicist’s inter- 
ests. Further sugges- 
tions will be most 
welcome. 


Analysis in 
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Flame ‘Temperature’ 


BY BERNARD LEWISt AND GUENTHER VON ELBEt 


Central Experiment Station, Bureau of Mines, and the Coal Research Laboratory, 
Carnegie Institute of Technology, Pittsburgh, Pennsylvania 


l. Introduction 


N this paper only flames produced incident to 
the use of fuels for heating and generation of 
power will be considered. If the fuel is in the 
gaseous state and is premixed with air or oxygen, 
after local ignition the burned and unburned 
gases are separated by a narrow reaction zone 
or flame front. If the fuel is discharged into the 
atmosphere combustion proceeds by interdiffu- 
sion of the gaseous or gasified fuel and the air. 
For organic fuels the latter flames are usually 
luminous because of formation of soot. This also 
applies to premixed flames with large excess 
of fuel. 
Other luminous phenomena to which the term 


flame has been applied, such as the cool flames 


observed in the slow oxidation of gaseous 
hydrocarbons and in the interdiffusion of alkali 
metals and halogens do not warrant discussion 
here. 


2. Theoretical and Experimental Flame 
Temperature 


The quantity temperature defines a state of 
statistical equilibrium and can be assigned only 
to a system in which such an equilibrium exists. 
The theoretical flame temperature thus corre- 
sponds to the following postulated state of 
complete statistical equilibrium in the burned 
gas: Collisions between the molecules and radia- 
tion, although continually changing the state 
of every individual molecule, maintain on the 
average the randomness of direction of motion 
and the constancy of (a) the percentage of 
molecules possessing a given velocity; (b) the 
percentage of molecules possessing a given 


* Published by permission of the Director, Bureau of 
Mines, U.S. Department of the Interior, and the Director, 
Coal Research Laboratory, Carnegie Institute of Tech- 
nology (not subject to copyright). 

+ Physical chemist, Explosives Division, Central Experi- 
ment Station, Bureau of Mines, Pittsburgh, Pa. 

t Coal Research Laboratory, Carnegie Institute of 
Technology, Pittsburgh, Pa. 
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quantum state of rotation, of vibration, or of 
electronic excitation, and (c) the concentrations 
of dissociation products. The theoretical flame 
temperature can be calculated' for a given 
content of total energy of the gas from heat 
capacities and dissociation equilibria obtained 
by statistical-mechanical methods from band 
spectroscopic and other data. It is reasonable to 
assume that in all zones of an actual flame the 
equilibrium (a) in the translational degrees of 
freedom is established with close approximation 
because of the rapidity of the collision sequences ; 
and since for perfect gases, which flame gases 
resemble, the temperature corresponding to this 
equilibrium is given by the gas law, an experi- 
mental flame temperature corresponding to the 
equilibrium in the translational degrees of 
freedom may always be obtained by determining 
the parameters of the gas law, pressure and 
molar volume. Whether the equilibrium (b) in 
the internal degrees of freedom and the chemical 
equilibrium (c) are also established with an 
approximation sufficient to eliminate measurable 
differences between theoretical and experimental 
flame temperatures remains to be investigated 
experimentally. 


3. Comparison Between Theoretical and 
Experimental Flame Temperatures in 
Explosion Experiments 

If an explosive gas mixture is ignited at the 
center of a spherical vessel, the resulting spherical 
flame does not come into contact with the cold 
walls of the vessel until combustion is complete, 
and the theoretical flame temperature can, in 
general, be calculated without difficult heat loss 
corrections. The maximum pressure established 
at the end of the explosion yields, on the other 
hand, the experimental flame temperature, the 
volume remaining constant. Experiments of this 
type have been carried out extensively for 


‘mixtures of hydrogen and oxygen with various 
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inert gases. It was found! that with the mon- 
atomic gases, argon and helium, as diluents the 
experimental and theoretical temperatures are 
identical within narrow limits of error, but that 
with diatomic gases, like oxygen or 
nitrogen, the experimental temperatures are 
appreciably higher than the theoretical. This in- 
dicates that the latter gases are slow to absorb 
energy in the internal degrees of freedom which 
they possess in distinction to the monatomic gases, 
and that equilibrium (b) is not established in 
the burned gas behind the flame front. The 
phenomenon has been referred to as excitation 
lag,! implying that the collisional excitation of 
the internal degrees of freedom, particularly 


excess 


those of vibration, is lagging behind the chemical 
reaction. Its occurrence suggests caution in the 
use of molecular infra-red radiation from flames 
for thermometric purposes, particularly if the 
emitting molecules are introduced into the 
mixture as diluent gas. 

Such excitation lag is also noticeable in experi- 
ments with sound waves.? A large rotational 
excitation lag has recently been reported in the 
rapid compression of diatomic gases at room 
temperature® and the rapid discharge of steam 
through an orifice.* 


4. The Line-Reversal Method of Determin- 
ing Flame Temperatures 


a. THEORY 


The line-reversal method of Féry® is based on 
the assumption that on introducing certain metal 
vapors into a flame statistical equilibrium is 
established between the electronic degrees of 
freedom of the metal atoms and the flame gases. 
The metal atoms thus emit and absorb their 
spectral lines as thermal radiators. 

The quantity measured by the line-reversal 
method is brightness, B, defined by the equation 


dE = BdAdtdw cos 6 (1) 


where dE is the energy radiated during the time 
dt from an element dA of the surface of a lumi- 
nous body within the solid angle dw at the 
angle @ to the normal. If the radiation is spec- 
trally unresolved B is the total brightness. If 
only the energy dE, contained in the portion dd 
of the spectrum (A being the wave-length) is 
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considered, then 
dE = B,dd\dAdtdw cos @. (2) 


By is called the spectral brightness. Fundamental 
optics as well as thermodynamics show that the 
brightness B or B, of a body remains constant 
in the path of any individual light ray emitted 
by the body regardless of the interposition of 
lenses or mirrors, provided that the optical 
system ideal one, that is, there is no 
diffraction or absorption by lenses and mirrors 
or reflection from lenses, and provided that the 
medium outside of the optical system is homo- 
geneous and nonabsorbing. 

The absorptivity, a, of any substance is defined 
as the fraction of incident radiant energy 
absorbed by it. The spectral absorptivity, ay, 
refers to a portion dd of the spectrum as above. 
A body that has an absorptivity of 1 is known 
as a blackbody. The emissivity, e, of any sub- 
stance is defined as the ratio of the brightness 
of the substance to the brightness of a blackbody 
at the same temperature. The spectral emissivity 
ey is defined correspondingly. 

For a thermal radiator of any kind and for 
any wave-length the relation between emissivity 
and absorptivity is given by Kirchhoff’s law: 

B,(nb) 


—=€,=Q), (3) 


By(b) 


is an 


where nb and 6 refer to nonblack and blackbodies, 
respectively. 

A blackbody obeys Lambert’s law, that is 
the brightness B or B, is independent of 6. It is 
a function of temperature alone, according to 
Planck’s radiation law: 

1 
— (4) 


By == CN 5 ’ 
ec? aD 1 


where ¢;>=Ch=1.40X10-" cal. cm? sec. (c is 
the velocity of light and h is Planck’s constant), 
and co=ch/k=1.433 cm deg. (k is Boltzmann’s 
constant), and T is expressed in °K. 

For short wave-lengths and throughout the 
visible up to about 0.74 and for temperatures up 
to about 3500°K, e%/*7>1 and Eq. (4) trans- 
forms to 

By) = :d-%e— 2 F, (5) 


known as Wien’s radiation equation. 
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The line-reversal method rests on the following 
principle: If a blackbody is placed behind a 
flame colored, for example, with sodium (which 
emits the yellow D-doublet, \=0.5890 —6yz), and 
a spectroscope is sighted through the flame onto 
the blackbody, there will be some temperature 
of the latter at which its brightness in the 
spectral region of the D lines equals the 
brightness of the light transmitted in this region 
through the flame plus the brightness of the D 
lines from the flame itself. At this temperature, 


Sp. D 


Le I Lj 
( Dyes pee 
ey, B C.R. 


4, 


Fic. 1. Diagrammatic sketch of line-reversal apparatus. 


therefore, only a continuous spectrum of the 
blackbody is visible in the spectroscope. At a 
lower temperature of the blackbody the sodium 
lines appear bright by contrast against the 
continuous spectrum. At higher temperatures 
this contrast reverses itself, the sodium lines 
appearing dark against the continuous spectrum. 
At the temperature at which this reversal takes 
place the following relation applies, if no light is 
reflected from the flame: 


Bow y t+Bow1—ap.y)) = Bo h)» (6) 


where the subscript D represents the wave- 
lengths of the D lines and f refers to the flame. 
Equation (6) can be transformed to 


Bory Bo 
=p (fs). (/) 


Bpw) 


This equation is of the form of Kirchhoff’s law 
(Eq. (3)); and if proof can be obtained of the 
thermal character of the sodium radiation the 
temperature of the flame is identical with the 
reversal temperature. 

The metal is usually introduced into the flame 
as a salt. The salt vaporizes and dissociates 
into metal atoms and other products. For the 
measurement of the experimental flame temper- 
ature as defined in Section 2 it is not necessary 
for the salt to vaporize completely or for dis- 
sociation equilibrium of the salt to be established. 
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It is only necessary that activation and deacti- 
vation of the metal atoms by collisions with 
other molecules be rapid enough to maintain 
the equilibrium concentrations of atoms in the 
electronic quantum states involved in the 
absorption and emission of the spectral line. 


b. EXPERIMENTAL TECHNIQUE 


Figure 1 shows a diagrammatic sketch of the 
line-reversal apparatus as used on stationary 
flamés. The image of a continuous light source 


‘(the comparison radiator, CR) is focused onto 


a portion of the flame above burner B by means 
of lens L,;. Diaphragm D defines the beam of light 
from the flame and comparison radiator which is 
brought to a focus on the slit of spectroscope S, 
by lens Le, the composite image completely cover- 
ing the slit. In practice it is not necessary to use 
a blackbody as comparison radiator; instead, a 
light source calibrated against a blackbody in the 
spectral range under consideration may be used. 
A convenient comparison radiator is a tungsten 
band lamp, the brightness of whose central 
portion can be determined as a function of the 
electric power input. The calibration can be 
made with an optical pyrometer which itself 
has been calibrated against a blackbody. The 
ordinary optical pyrometer utilizes only red 
brightness by use of a red filter, \=0.665y, but 
the calibration can be extended to other parts 
of the tungsten spectrum by using spectral 
emissivity data® and Eqs. (3) and (5). Other 
comparison radiators that have been used are 
a carbon filament lamp,® a Nernst lamp,” ® and 
for very high temperatures a tungsten arc® and 
carbon arc!’ with a rotating sector or absorption 
glass for diminishing the brightness. If lens 
L,; is used between the comparison radiator and 
the flame, a correction for absorption and 
reflection due to the lens must be made." 

The metal may be introduced into the flame™ 
either as a spray of a solution of the salt or as a 
dry dust. The latter may be introduced directly 
into the flowing gas by mechanical means or 
by passage of a spark or arc between salt 
electrodes containing a metal core™ or between 
electrodes immersed in powdered salt. In some 
investigations a piece of salt was placed on the 
grid of the burner and was vaporized by the 
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heat of the flame. The latter method permits 
coloring of a narrow portion of the flame. 
It can be accomplished by the spray or dust 
method if the salt is introduced with a portion 
of the gas stream passing through a separate 
smaller tube placed within the burner. 

Modifications of the technique of observing 
the reversal point have been proposed™ utilizing 
the contrast effects afforded by uneven bright- 
ness along the observed portion of the tungsten 
strip. 


c. VALIDITY OF THE METHOD 


There is no instance in the literature in which 
the gas law parameters have been determined 
simultaneously with the line-reversal tempera- 
ture. Therefore, within the definition of experi- 
mental flame temperature (Section 2) the line- 
reversal method has never been tested. With 
Bunsen burners this could be done in principle 
by measuring the refraction of the gas stream 
on passing through the boundary of the inner 
cone using the relationship’ 


p, tang 
= - —, (8) 
p, tana 


where p, and p are the densities of the unburned 
and burned gases, respectively, a the angle at 
which the gas stream enters and £ the angle at 
which it leaves the reaction zone. Accurate 
measurements of 8 are very difficult. 

There is, however, a wealth of other material 
which has established the validity of the line- 
reversal method. Using a large Méker burner and 
coal gas-air mixtures, Bauer’ found identical 
reversal temperatures with a number of different 
salts of sodium and salts of other metals with 
characteristic line emissions in all parts of the 
visible spectrum. It does not seem possible to 
interpret this result in any other way than that 
statistical equilibrium is established in all cases, 
either completely or at least between the elec- 
tronic degrees of freedom of the metal atoms and 
the translational degrees of freedom of the gas. 
Since the mixture compositions were not known 
well enough to permit comparison between the 
theoretical calculated temperatures and the re- 
versal temperatures, no estimate can be made of 
the extent to which equilibrium was attained in 
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these experiments. Kohn*® confirmed Bauer’s 
result of the reversal temperature on the nature 
of the salt, using again a Méker burner and coal 
gas-air mixtures of unknown composition. Fur- 
thermore, Kohn compared the reversal tempera- 
ture with the temperature of a wire immersed in 
the flame and compensated electrically for radi- 
ation using a procedure described by Schmidt.'® 
A wire, inserted in a flame, emits radiation in a 
continuous spectrum for which the flame is 
essentially transparent ; it reaches a temperature, 
lower than that of the gas, at which the flow of 
energy from the gas to the wire equals the 
energy radiated from the wire. The wire is now 
heated by an electric current to a temperature at 
which there is no further flow of energy from the 
gas to the wire. The exact amount of current can 
be determined from knowledge of the radiated 
energy per unit length of wire as a function of the 
wire temperature and increasing the wire temper- 
ature in the flame until the electric power input 
equals the radiant energy. The temperature so 
determined was found to be identical with the 
line-reversal temperature measured in the neigh- 
borhood of the wire. Thus, the same equilibrium 
exists between the gas and the solid metal as 
between the gas and the vaporized metal atoms. 

Calculation of theoretical flame temperatures 
for known mixtures is complicated by heat 
losses to the burner grid and also to the sur- 
roundings by radiation and intermixing. In the 
case of excess fuel, secondary combustion occurs 
in the outer parts of the flame. There exists, 
however, a zone of considerable height and width 
above the grid in which the temperature is quite 
uniform and in which, therefore, the disturbance 
due to contact with the atmosphere is negligible, 
and only losses to the grid and by radiation are 
important. Reversal temperatures in this zone of 
Méker flames of known coal gas-air composition 
and varying flow velocities of the gas were 
determined by Minkowski, Miiller and Weber- 
Schafer.'® Comparison of their values with the 
theoretical flame temperatures calculated for 
adiabatic combustion shows that for sufficiently 
high flow velocity the reversal temperature ex- 
ceeds the theoretical temperature (see Fig. 2).* In 


* Two limiting curves for dry and saturated gases were 
calculated since the water vapor content was not known. 
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Stor. 











view of Bauer’s and Kohn’s results establishing 
the validity of the line-reversal method, this 
indicates the existence of excitation lag in the 
flame gas similar to that found in explosion 
experiments described in Section 3. Figure 2 also 
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Fic. 2. Comparison of theoretical flame temperatures of 
coal gas-air mixtures with experimental flame temperatures 
determined by the line-reversal method. — - — theoretical 
curve, dry mixtures. ———— theoretical curve, moist 
mixtures. @ sodium line-reversal flame temperatures with 
gas velocity in centimeters per second (Minkowski, et al.). 


shows that the gas velocity required for the 
reversal temperature to transcend the theoretical 
temperature becomes larger toward stoichio- 
metric mixtures. This is believed to be in part a 
consequence of the higher burning velocity of 
mixtures nearer the stoichiometric because the 
higher the burning velocity the closer the flame 
burns to the grid for a given flow velocity and the 
greater the heat loss. Part of the effect must also 
be attributed to the higher flame temperatures of 
such mixtures causing more rapid heat transfer. 

Reversal temperature measurements by 
Kaveler and Lewis" on slow burning natural gas- 
air flames agree very closely with the theoretical 
temperatures for adiabatic combustion (Fig. 3). 
Only in the neighborhood of the stoichiometric 
mixture do the reversal temperatures fall some- 
what below the theoretical curve which may be 
explained in the same way, namely, greater heat 
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loss to the grid for higher burning velocities. If 
oxygen is substituted for the atmospheric nitro- 
gen in lean mixtures the reversal temperatures 
are, on the average, somewhat above the theo- 
retical curve (Fig. 4) which is another indication 
of excitation lag, it being known from explosion 
experiments that this effect is considerably larger 
for oxygen than for nitrogen. 

In all the investigations mentioned the experi- 
menters were careful to color only the portion of 
the flame being investigated. If the entire flame is 
colored the reversal temperature should lie be- 
tween values for the gas at the center and the gas 
at the boundary. The latter is cooler than the 
former in lean mixtures, and hotter, owing to 
secondary combustion, in rich mixtures. Some 
early work of Kurlbaum and Schulze” which 
seemed to indicate a dependence on the nature of 
the salt was shown by Bauer’ to be in error be- 
cause of the use of totally colored flames. 

In flames of organic fuels the reaction zone is 
colored blue or green (Swann bands) depending 
on the mixture composition. This emission is due 
to C—H or C—C (Swann bands) radicals, 
respectively. Since these radicals are formed as 
intermediates in the combustion reaction and 
probably are excited in the act of formation, it 
might be expected that their electronic quantum 
states are not in statistical equilibrium. This is 
confirmed by observations of Bauer’ who applied 
a modified line reversal method to these bands 
and obtained temperatures, obviously fictitious, 
which were many hundreds of degrees higher 
than the possible maximum temperature of the 
flame. This result has been confirmed inde- 
pendently by Tawde and Patel.'® 


5. Determination of Flame Temperature by 
Wires and Thermocouples 


If a wire or thermocouple is inserted in a flame 
for thermometric purposes allowance must be 
made for the temperature difference between 
wire and gas due to radiation (see Section 4c). 
This may be done in two ways: (1) By using a 
series of progressively thinner wires and ex- 
trapolating the measured wire temperatures to 
zero thickness. (2) By compensating the radiation 
by means of electrical heating. The latter method 
requires determination of the energy radiated 
from the wire per unit length and time as a 
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function of its temperature. This may be done 
either by measuring the radiation with a thermo- 
pile for a series of pyrometrically determined 
wire temperatures’ (these measurements may 
be performed with the wire in any transparent 


medium, such as the open atmosphere*), or by 


measuring the electric power input for a series of 
pyrometrically determined temperatures with 
the wire placed in a high vacuum so that con- 
duction and convection losses are eliminated.® 

In the various suction pyrometers!® used in 
furnace work the thermocouple is enclosed in a 
tube through which the furnace gas is passed rap- 
idly enough to bring the walls of the tube to the 
temperature of the gas. There is thus no radiation 
loss from the thermocouple to the colder furnace 
walls. Sometimes the walls of the tube are heated 
electrically with the current so adjusted that the 
thermocouple reading remains the same whether 
the gas is flowing or not. 


6. Determination of Flame Temperature by 
the Measurement of Brightness and 
Absorptivity in the Infra-Red 

The temperature of a thermal radiator can 
be determined by direct measurement of its 
brightness and absorptivity in some spectral 


2000 





Stoichiometric ~ 


~ 


Flame temperature, °C 


1600 





6 8 10 11 12 
Pittsburgh natural gas, percent in air 


Fic. 3. Highest flame temperatures of Pittsburgh natural 
gas-air mixtures observed in center of Méker burner flame 
by sodium line-reversal method. — theoretical flame 
temperatures. © observed flame temperatures (Kaveler 
and Lewis). 





range. The relation between these quantities is 
given by Eqs. (3) and (4). Such measurements 
were made by Schmidt" for wave-lengths be- 
tween 2.74 and 4.4u (bands of HeO and COs) 


* Convection and conduction losses do not affect these 
measurements as has sometimes been erroneously supposed. 
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Fic. 4. Highest flame temperatures of Pittsburgh 
natural gas-oxygen mixtures observed in center of Méker 
burner flame by sodium line-reversal method. ———— theo- 


retical flame temperatures. O observed flame temperatures 
(Kaveler and Lewis). 


and by Bauer’ for wave-lengths in the region of 
25u (residual rays of fluorite). Both investigators 
used large Méker burners and coal gas-air mix- 
tures. Schmidt also measured the temperature of 
the flame by a wire in the manner described 
above and found fairly good agreement between 
the temperatures obtained by both methods. 
Bauer found similar agreement between his 
measured infra-red temperature and the line- 
reversal temperature. It is not easy to appreciate 
the significance of these experiments because 
infra-red absorption and emission take place 
throughout the entire depth of the flame on which 
the instrument is sighted and it is not possible to 
single out a narrow zone as in the line-reversal 
and wire methods. Thus, the temperature ob- 
tained by the infra-red method should be inter- 
mediate between extreme values existing in the 
flame, as would be the case with the line-reversal 
method if totally colored flames were used. 


7. Measurement of Temperature and Emis- 
sivity of Soot-Forming Luminous 
Flames 


The slowness of the combustion process in 
soot-forming, luminous flames of the pre-mixed 
type with rich mixtures or of the diffusion type 
insures the absence of any appreciable excitation 
lag effect, particularly in view of the experience 
of Kaveler and Lewis'* with natural gas-air mix- 
tures burning relatively much faster. Thus, the 
luminous flame has a temperature correspond- 
ing to the definition of the theoretical flame 
temperature—it is a thermal radiator, the bulk 
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of the radiation being emitted by the soot par- 
ticles as a continuous spectrum. The soot par- 
ticles, being cooled by radiation more strongly 
than the surrounding gas, are at a somewhat 
lower temperature than the latter. On the basis 
of data on particle size and heat transfer of solid 
bodies in gaseous media, Schack*® has estimated 
that the temperature difference between the soot 
particles and the gas should be substantially less 
than 1°C at usual flame temperatures. 

For measuring the temperature of a luminous 
but not opaque flame Kurlbaum* proposed a 
method similar to the line-reversal method. An 
optical pyrometer is sighted on a blackbody with 
and without the flame between them. The tem- 
perature of the blackbody is adjusted until 
its brightness, as seen through the pyrometer, 
is the same as that observed with the flame 
interposed. Equation (6) applies which means 
that both flame and blackbody are at the same 
temperature, given by the pyrometer reading. 
Other comparison radiators may be substituted 
for the blackbody. This method is not easily 
applicable to industrial furnaces. 

Since the presence of soot particles greatly 
increases the heat transfer by radiation in 
furnaces much thought has been given to the 
determination of emissivity of luminous flames 
from temperature measurements. Schack,”° after 
reviewing the data on light absorption of soots 
from various sources, came to the conclusion 
that the spectral absorptivity of the flame can be 
represented by the absorption law for homo- 
geneous absorbing media, namely 


a,=1—e*r" (=e, see Eq. (3)) (9) 


where L is the thickness of the absorbing layer 
and k, is the absorption coefficient which, over a 
limited range of wave-lengths can be represented 
by the equation 
k 
k=—, (10) 
A 
where k and a are constants. Introducing the 
spectral brightness temperature T, of the flame, 
which is the temperature of a blackbody of equal 


spectral brightness, one obtains from Eqs. (3) 


and (5) 


———=— Ing, (11) 
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where 7 is the true temperature of the flame. If 
7, is measured for a given wave-length and T is 
determined independently, then from Eqs. (9) 
and (11) the absorption strength k,L may be 
calculated ; and since values of a are known over 
the entire spectral range, e, is also known as a 
function of the wave-length. The desired quantity, 
the total emissivity of the flame, can then be 
calculated from  Kirchhoff’s law for total 
brightness, 


S0%exBynydd 


e= —, (12) 
So” Bywydd 
By») being given as a function of \ by the Planck 
equation (4). 
Schack suggested measuring the brightness 
with an optical pyrometer and the true tempera- 
ture by some other method such as a suction 
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Fic. 5. True temperature and absorption strength of 
luminous flames as functions of red and green brightness 
temperatures (Hottel and Broughton). 


thermometer. But, as he himself pointed out, the 
latter method does not allow a precise determi- 
nation of the usually small difference between 7) 
and 7, mainly because of the use of two entirely 
different types of instruments. Hottel and 
Broughton” overcame this difficulty by em- 
ploying the optical pyrometer exclusively and 
determining 7) in two spectral regions, the red 
and green. Two equations of the type of 11 may 
now be written with 7), = 7 peqg and Treen. Hence, 
T can be eliminated, otherwise the procedure for 
calculating e remains as before. Hottel and 
Broughton made further studies of the best 
values of a to be used. They are a=1.39 up to 
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A=0.84 and a=0.95 from A=0.8u to 10n. 
To eliminate the tedious calculations Hottel 
and Broughton constructed two graphs (Figs. 
5 and 6) from which the true temperature and 
total emissivity may be read knowing Tz and T¢. 
From the value of Tr on the abscissa of Fig. 5 
one goes vertically to the proper value of 
(T'g—Tr), then reads the true temperature T by 
interpolation on the semivertical lines and kL 
on the ordinate. With this value of RL, e is read 
from Fig. 6 for any true flame temperature. 

Another method adapted to industrial work, 
namely color pyrometry,”* makes use of the fact 
that for many nonblack radiators (so-called gray 
bodies) the ratio of the spectral brightness for two 
wave-lengths is about the same as fora blackbody. 
In its simplest form the color pyrometer consists 
of a bichromatic wedge-shaped light filter that 
transmits only red and green. The coefficients of 
extinction of the two colors differ; as the thick- 
ness of the wedge is increased green is absorbed 
more completely than red. When a radiating 
body is observed through this type of filter it 
appears green when viewed through the thin end 
and red through the thick end. At a certain point 
in the wedge these two colors, being comple- 
mentary, merge into white. The position of this 
point in the wedge depends on the ratio of the 
red and green brightnesses of the radiating body, 
i.e., on its temperature, and the length of the 
wedge can thus be calibrated in terms of temper- 
ature. Other more elaborate instruments have 
been designed with compensating devices for 
nongray bodies and for overcoming physiological 
errors. 


8. Concluding Remarks 


Experimental flame temperatures significantly 
higher than theoretical temperatures have been 
observed only in fast flames such as hydrogen 
and coal gas flames. In furnace flames there 
appears to be no reason to expect appreciable 
deviations from statistical equilibrium except 
perhaps near the orifice in pre-mixed flames— 
disturbances such as heat transfer are too slow 
to affect the statistical equilibrium. Thus, there 
is no anomaly to be dealt with and temperature 
measurements should give the true flame temper- 
ature that corresponds to the theoretical flame 
temperature as defined in Section 2. 
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Fic. 6. Total emissivity of luminous flame from its 
true temperature and absorption strength (Hottel and 
Broughton). 


In internal combustion engines where the flame 
is greatly speeded up by turbulence, appreciable 
excitation. lag may occur, but owing to the com- 
plexity of the process there seems to be no 
practical way at present to determine its 
contribution to the pressures developed. Rass- 
weiler and Withrow* determined a temperature 
of 2060°C at the firing end of an Otto cycle 
engine using the line reversal method. This 
temperature was calculated later by Hottel and 
Eberhardt® to be 2071°C on the basis of adia- 
batic combustion. The measured temperature is 
undoubtedly lowered by cooler gas layers near 
the observation window. 

Flame temperatures of various combustibles 
and the change of flame temperature with the 
air-fuel ratio can be compared from calculations 
on the basis of adiabatic combustion. These 
calculations are tedious and have been done with 
accurate data only for a few cases.' In most in- 
stances where less accurate older data have been 
used the results are in error by not more than 50° 
and may still be used as a basis of comparison of 
various mixtures. Lacking calculated values, the 
large amount of data of line-reversal tempera- 
tures may prove useful although these data were, 
as a rule, obtained with totally colored flames. 
Such line-reversal temperatures are uniformly 
lower than the theoretical temperatures by 
amounts ranging from 30 to 80°C. Extensive 
tables have been compiled by Lewis and von 
Elbe.*® 

In gases where statistical equilibrium is 
very strongly disturbed, such as in discharge 
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tubes and in the hydrogen torch, the term 
temperature loses much of its practical signifi- 
cance although it may still be a definable 
quantity on the basis of the translational energy 
of the gas. Pirani and Rompe” have proposed 
to denote the degree of freedom or the molecular 
state that primarily acquires the energy added 
to the system as “‘energy flow trap.”’ For ordinary 
flame gases this term does not seem to be 
usefully applicable. 


9. Summary 


The discussion is confined to flames produced 
incident to the use of fuels for heating and 
generation of power. The theoretical flame 
temperature corresponding to complete statisti- 
cal equilibrium in the burned gas is distinguished 
from the experimental flame temperature meas- 
ured by the gas law equation and corresponding 
to equilibrium in the translational degrees of 
freedom. Experiments on fast flames in closed 
vessels show that the latter may be higher than 
the former due to excitation lag in the internal, 
presumably vibrational, degrees of freedom. 





The theory and technique of the line-reversal 
method are described and the experimental evi- 
dence is interpreted as indicating that the method 
yields the experimental flame temperature as 
defined above. The technique of measuring 
flame temperatures by wires and thermocouples 
is described. In measurements of flame tempera- 
tures by infra-red radiation temperature varia- 
tion within the flame is a disturbing factor. 

Soot-forming luminous flames may be regarded 
to be in complete statistical equilibrium, the 
temperature difference between soot particles 
and the surrounding gas being negligible. 
Methods of determining the temperature and 
emissivity of such flames are described. In 
industrial furnaces the differences between the 
experimental and theoretical temperatures should 
be negligible and there is thus no reason to 
question the existence of a true temperature. 
In internal combustion engines excitation lag is 
likely to exist but no experimental confirmation 
is available. The term “energy flow trap’’ does 
not seem to have useful application in technical 
combustion processes. 


Bibliography 


1. B. Lewis and G. von Elbe, Combustion, Flames and 
Explosions of Gases (Cambridge University Press, 1938). 

2. See, for example, H. O. Kneser and V. O. Knudson, 
Ann. d. Physik 21, 682 (1935), and other references 
mentioned there. 

3. B. Lewis and G. von Elbe, J. Chem. Phys. 7, 197 
(1939). 

4. B. V. Korvin-Kroukovsky, J. Frank. Inst. 221, 99 
(1939). 

5. Ch. Féry, Comptes rendus 137, 909 (1903). 

6. W. E. Forsythe and A. G. Worthing, Astrophys. J. 
61, 140 (1925). 

7. E. Bauer, Recherches sur le rayonnement, Théses 
(Gauthier-Villars, Paris, 1912). 

8. H. Kohn, Ann. d. Physik 44, 749 (1914). 

9. E. Griffiths and J. H. Awbery, Proc. Roy. Soc. A123, 
401 (1929). 

10. F. Henning and C. Tingwaldt, Zeits. f. Physik 48, 805 


(1928), H. H. Lurie and G. W. Sherman, Ind. Eng. Chem. 


25, 404 (1933). 

11. For details see G. W. Jones, B. Lewis, J. B. Friauf and 
G. St. J. Perrott, J. Am. Chem. Soc. 53, 869 (1931). 

12. Cf. G. Ribaud, Y. Laure and H. Gandry, Inst. of 
Fuel, Symposium on Gas Temperature Measurement, Vol. 
12, March, 1939, pp. 18-30. 

13. H. H. Kaveler and B. Lewis, Chem. Rev. 21, 421 
(1937). 


706 


14. Gouy, Ann. Chim. Phys. (5) 18, 27 (1879); cf. refer- 
ence 1. 

15. H. Schmidt, Ann. d. Physik 29, 1027 (1909). 

16. R. Minkowski, H. G. Miiller and M. Weber-Schiafer, 
Zeits. f. Physik 94, 145 (1935). 

17. F. Kurlbaum and G. Schulze, Verh. d. D. Phys. Ges. 
5, 427 (1903). 

18. N. R. Tawde and J. M. Patel, J. Univ. of Bombay 
6, 29 (1937). 

19. Cf. Inst. of Fuel, Symposium on Gas Temperature 
Measurement, Vol. 12, March, 1939. 

20. A. Schack, Zeits. f. tech. Physik 6, 530 (1925). 

21. F. Kurlbaum, Physik. Zeits. 3, 187 (1902). 

22. H. C. Hottel, and F. P. Broughton, Ind. Eng. Chem. 
(Anal. Ed.) 4, 166 (1932). 

23. G. Naeser, Inst. of Fuel, Symposium on Gas Tem- 
perature Measurement, Vol. 12, March, 1939, p. 38; W. 
Bowen, tbid., p. 75. 

24. G. Rassweiler and L. Withrow, J. Soc. Aut. Eng. 36, 
125 (1935). 

25. H.C. Hotteland J. E. Eberhardt, Chem. Rev. 21, 439 
(1937). 

26. Reference 1, pp. 398-401. 

27. M. Pirani and R. Rompe, Trans. Electrochem. Soc. 
69, 417 (1936); cf. J. G. Bennett and M. Pirani, Inst. of 
Fuel, Symposium on Gas Temperature Measurement, Vol. 
12, March, 1939, p. 1. 


JOURNAL OF APPLIED PHYSICS 























Resumeés of Recent Research 








Thermal Conductivity 
of Metals 


From the standpoint 
of modern theory, the 
expression for thermal 
conductivity of a metal should contain two 
additive terms, one due to the lattice and one 
due to the electron conductivity. A. H. Wilson! 
gives the equation as 

1 nr? ntkyT 
k=-—pClu+--—- —., (1) 
3 3 m 
k is the thermal conductivity; p, the density; C, 
the atomic heat; /, the electron mean free path 
in the lattice; u, the velocity of elastic waves 
(sound) in the lattice; ”, the number of free 
electrons per unit volume; 7, the “relaxation 
time’ (average time between collisions); Ro, the 
molecular gas constant; m, the mass of the elec- 
tron; 7, the absolute temperature. 
If we write the equation as 


k 1 T" ntkyT 
-=—lu+ (2) 
pC 3 3 mpC 
it is in the form 
k 1 
=K—+K’. (3) 
eC T 


This assumes that 7 « 7~! and lu « 7—' and that p 
and C compensate as to their temperature effects. 

The equation in this latter form is experi- 
mentally verified by data on tin, lead and zinc.* 


Values of k/ pC, plotted as ordinates against 1/7 


as abscissa, give straight lines with intercepts on 
the vy axis and these intercepts turn out to be the 
values of k/ pC for the molten state. 

The conductivity for the molten state includes 
the residual atomic conductivity as well as that 
due to electrons. The conductivity of an insulat- 
ing liquid is very small (1 percent, or less, of 
that of a conductor). Hence the last term on the 
right is, very likely, almost entirely electron con- 
ductivity. The temperature effect for the solid 
is seen to be entirely, or nearly entirely, due to 


' A. H. Wilson, Semiconductors and Metals (Macmillan, 
New York, 1939). 

* This note is a resume of papers by C. C. Bidwell in 
Phys. Rev. 56, 594-598 (1939); 58, 561-564 (1940). A 
precise method is there described for making thermal con- 
ductivity measurements at high temperatures, applicable 
to either solid or molten metals. 
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the elastic lattice. The data indicate little or no 
temperature effect in the liquid state. 

Thermal conductivity measurements on other 
metals offer much greater difficulty than was the 
case for those here reported, which are fairly low 
melting point metals. However, from the law 
indicated, one may venture to predict the ther- 
mal conductivity for the molten state for any 
metal provided enough data are available on the 
solid state to establish the straight line and 
obtain the intercept. 

Measurements were not carried below the tem- 
perature of liquid nitrogen (—195°C) and it is 
not expected that these relations should hold at 
extremely low temperatures. 


Young’s Modulus of 
Beta-Brass Crystals 


The marked changes 
in the physical proper- 
ties of certain ordered 
alloys as they go through the order-disorder 
transition have been measured extensively for 
polycrystalline material, but little work has been 
done with single crystals. Recently,' however, the 
change in Young’s modulus of seven single 
crystals of beta-brass, an approximately equi- 
mixture of copper and zinc has been measured as 
a function of temperature. The method used was 
that of the composite piezoelectric oscillator. 
From such measurements it is possible to obtain 
Young’s modulus for any direction in the crystal 
lattice. Plots of the reciprocal of Young’s modulus 
in the three principal directions, [100], [110], 
and [111], (cube edge, face-diagonal, and body- 
diagonal, respectively), show that the effect of 
disorder is mainly felt in the first two directions, 
while in the third direction the elasticity shows an 
approximately linear increase with increasing 
temperature. This is very surprising in view of 
the fact that in almost every other known 
material just the reverse is true. In the [111 ] 
direction the behavior of Young’s modulus is like 
that for an isotropic solid, while in the [110] 
direction there is almost zero change in 1/E. An 
unusual degree of anisotropy in Young’s modulus 
at room temperature, as shown by the large 
ratio, 8.2, of maximum to minimum modulus in 
the [111] and [100] directions, was noted, and 
this may account for the peculiar temperature 
dependence of Young’s modulus in different 
directions in the crystal. 


1 J. Rinehart, Phys. Rev. 58, 365 (1940). 


707 











Here and There 





Increased Facilities for Optical Work 


In order to provide facilities for the production of optical 
instruments and optical parts required by the National 
Defense Program, the Spencer Lens Company, Scientific 
Instrument Division of American Optical Company, has 
let a contract to The Gillmore-Carmichael-Olson Company 
for the construction of additional buildings at the Eggert 
Road site, Buffalo, New York, Mr. B. H. Witherspoon, 
president of the Company, announces. One unit of ap- 
proximately the same size as the present Mechanical Parts 
Plant, completed last summer, will house certain prism 
and lens manufacture. A connecting unit two and three 
stories high will be utilized for the assembly of special 
contract work. A third unit will provide expansion for the 
mechanical parts operations. It is planned to have the new 
unit in operation not later than February 1, 1941. 


* 


National Defense Products 


General Electric leads American industry in variety of 
products being made for the national defense program, ac- 
cording to a statement made by Philip D. Reed, General 
Electric Board Chairman, in a recent address before the 
Chamber of Commerce and the Junior Chamber of Com- 
merce of San Francisco. “‘A significant portion of the 
General Electric Company's productive capacity will be 
fully occupied with rearmament work,” Mr. Reed pointed 
out. He also stated, however, that he believed the Com- 
pany’s normal output of useful goods can be maintained 
and in certain departments substantially increased. The 
continuation on the part of industry of the production of 
useful goods is important, in Mr. Reed's opinion, if inflation 
is to be prevented and American standards of living 
preserved. 


* 


National Chemical Exposition 


New discoveries and developments in the service of 
chemistry and its application and progress in many fields 
of activity will be discussed by foremost authorities in 
America at an Industrial Chemical Conference, which will 
be a feature of the National Chemical Exposition to be held 
December 11 through 15 at the Stevens Hotel in Chicago. 

“We have arranged for this Industrial Chemical Con- 
ference an excellent program of papers in the nature of 
symposia on ‘Newer Developments in Chemistry and 
Chemical Engineering,’” said Dr. R. C. Newton, chairman 
of the show committee. ‘The program is timed to afford 
those who wish to hear the talks in which they are inter- 
ested ample opportunity also to view the show.” Technical 
sessions for this conference have been arranged for Thurs- 
day afternoon and evening, December 12; Friday forenoon, 
December 13, and Saturday afternoon, December 14. 
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New Calcite Source Discovered 


A new source of calcite, or Iceland spar, a valuable 
optical mineral heretofore imported from abroad, was re- 
vealed recently when officials of the Bausch & Lomb 
Optical Company announced they had contracted for the 
output of a new mine in the San Pedro mountains, north- 
west of Santa Fe, New Mexico. The discovery by a Mexican 
prospector of the new source of supply and the interest of 
E. M. Stanton, a Chicago optometrist, led to a subsidy by 
Bausch & Lomb through which operations have been 
carried on. 

The scarcity of optical calcite has caused a world-wide 
search for many years, since the mineral is essential in the 
construction of all polarizing instruments. With the 
flooding of the Iceland mine at Helgustadir during the 
World War, this source which formerly supplied the world 
has been unproductive of fine crystals. Meager supplies 
have come from various countries, the most recent from 
Spain and South Africa, but the total volume offered for 
sale has been insignificant for many years. 

Officials reported that the new mine had delivered more 
than five hundred pounds of fine calcite crystals within a 
period of three months, many of them weighing as high as 
forty pounds. An analysis of the output of the new mine 
indicates that the United States now has a supply of calcite 
of high optical quality which will serve the country for 
years. 


* 


Company Expansion Planned 


The Westinghouse Electric and Manufacturing Company 
will spend more than $8,500,000 during the year ending 
next June, to increase its emergency and peacetime pro- 
duction facilities. The Company’s expansion authorization 
provides some $5,000,000 for machine tools and equipment, 
in addition to the construction of 15 new buildings in six 
different states. Another $6,000,000 has been authorized for 
maintenance of the expanded plant facilities. 

Since last June, Westinghouse has completed five build- 
ings; nine are under construction, and another is con- 
templated. Two structures now nearing completion are a 
temporary building for production of radio equipment for 
the Government at the Westinghouse Radio Division in 
Baltimore, Maryland, and a service plant at Bridgeport, 
Connecticut. Two additional buildings for Government 
radio production and a service plant are under construction 
at Baltimore. Three factory additions and a warehouse are 
being built at the Westinghouse Merchandising Division, 
Mansfield, Ohio. 

A $500,000 office building at the Westinghouse Steam 
Division, South Philadelphia, and additional factory space, 
including a new printing plant, at the Westinghouse 
Micarta Plant, Trafford, Pennsylvania, are under way. A 
$300,000 district headquarters and service building has 
been authorized for Atlanta, Georgia. 
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University of Maine Faculty Promotions 


The following promotions in the Department of Physics 
at the University of Maine have been made this year: 

Professor Clarence E. Bennett, formerly Associate Pro- 
fessor and Department Head, to Professor and Department 
Head; 

Dr. Arthur O. Williams, Jr., from instructor to Assistant 
Professor; 

Mr. Frederick B. Oleson from assistant to part-time 
instructor. 


* 


Meteorology 


Commander F. W. Reichelderfer, of the U. S. Weather 
Bureau, has made a statement to the effect that only three 
institutions in the country now have an advanced meteoro- 
logical course, and these have each graduated from ten to 
twenty men a year. Most of the men have been Army and 
Navy officers sent there by the military service. There have 
been only about ten to twenty men each year for civilian 
occupations. This is due to the fact that the opportunity 
for employment has been limited. However, he states that 
the new five-day weather forecasting system, initiated 
recently by the Weather Bureau, is creating a demand for 
meteorologists.—Science 


* 
Mathematical Biophysics 


The table of contents of the December issue of the 
Bulletin of Mathematical Biophysics, published by the 
University of Chicago Press, includes the following articles: 

Some General Theorems on the Motion of Incompressible Viscous 
Fluids, by G. YouNG 

A Neural Mechanism for Discrimination. III. Visually Perceived 
Lengths and Distances, by A. S. HOUSEHOLDER 


Deformation of Shell Shaped Cellular Aggregates: Application to 
Gastrulation, by N. RASHEVSKY 


Some Considerations on Mathematical Molecular Biophysics, by N. 
RASHEVSKY 

Communications for the editor and manuscripts should 
be addressed to N. Rashevsky, Editorial Office of the 
Bulletin of Mathematical Biophysics, 5822 Drexel Avenue, 
Chicago, Illinois. 


* 


Color Television 


From the Columbia Broadcasting System comes an 
announcement of its laboratory success in developing color 
television. It is hoped that the manufacture of commercial 
equipment in this field will be started within the next few 
months, although, since Columbia is not itself a commercial 
manufacturer of television apparatus, it is not possible for 
the Company to say how long it will take for the actual 
marketing of the necessary devices to be begun. 


* 


Dr. Hans Bethe, Professor of Physics at Cornell Uni- 
versity, has been granted a sabbatical leave from February, 
1941, until February, 1942. He will spend the spring 
semester of 1941 at Columbia University, where he has 
been invited to lecture.—Science. . 
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35,000-Foot Parachute Jump 


For the purpose of solving some of the problems of 
stratosphere flying, a six-and-a-half-mile parachute jump, 
to be made in the interest of science by amateur Arthur H. 
Starnes, was being planned when this material went to 
press. In the laboratories of the Armour Research Founda- 
tion of Illinois Institute of Technology, scientists of the 
University of Chicago and Northwestern’ University 
gathered on September 23 to conduct final tests pre- 
liminary to the proposed 35,000-foot leap to be undertaken 
sometime within the two weeks following that date. 
Unique in the history of science and aviation, the tests 
simulated conditions of temperature and wind velocities 
which Starnes was expected to encounter as he hurtled 
earthward at 200 miles per hour. Fully clothed in jumping 
regalia, replete with 103 pounds of scientific instruments, 
Starnes was put through the cold wind tunnel tests. 

In addition to the interest of science in the jump planned, 
it was revealed that much of the data and information 
obtained thus far through the researches carried on by 
experts from Chicago and Northwestern Universities will 
be turned over to the officers of the Army and Navy, whose 
interest in them is based on the increased importance of 
parachute troops. 


* 


Executives Receive Promotions 


W. Richison Schofield, formerly chief engineer of Leeds 
& Northrup Company, Philadelphia, has recently been 
appointed Director of Engineering. John W. Harsch, 
Assistant Chief Engineer, has been advanced to Chief 
Engineer, and Mr. Harsch’s former position has been filled 
by the promotion of John F. Quereau. 


* 


Calendar of Meetings 
October 
28-Nov. 2 Conference on Applied Nuclear Physics, Cambridge, 
Massachusetts 
November 


11-15 American Petroleum Institute, Chicago, Illinois 
15-16 Acoustical Society of America, Chicago, Illinois 
22-23 American Physical Society, Chicago, Illinois 


December 
2—- 4 American Institute of Chemical Engineers, New Orleans, 
Louisiana 


2-7 Highway Research Board of National Research Council, 
Washington, D. C. ; i ; 

3-6 American Society of Mechanical Engineers, New York, New 
York 

20-21 American Physical Society, Pasadena, California ; 

26-28 American Physical Society, Philadelphia, Pennsylvania 

26-28 Geological Society of America, Austin, Texas ; ' 

28-30 American Association of Physics Teachers, Philadelphia, 
Pennsylvania ; “7 

27-Jan. 2 American Association for the Advancement of Science, 
Philadelphia, Pennsylvania 


February APs. 
21-22 Joint meeting, American Physical Society and Optical Society 
of America, Cambridge, Massachusetts 
April 
7-11 American Chemical Society, St. Louis, Missouri 
May ; 
1- 3. American Physical Society, Washington, D. C. 
June 


20-21 American Physical Society, Providence, Rhode Island 
23-27. American Association for the Advancement of Science, Durham, 
New Hampshire 
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Biophysics Fellowship 


Union College, Schenectady, New York, made recent 
announcement of a fellowship in biophysics for 1940-41. 
Further information may be obtained from President Dixon 
Rayan Fox, Union College, Schenectady, New York. 


Conference on Electrical Insulation 


The Conference on Electrical Insulation of the National 
Research Council (Division of Engineering and Industrial 
Research) is holding its annual meeting in Washington, 
DD. C., on October 31 and November 1 and 2. 


Contributors to This Issue 








Engineering. 
C. S. Barrett ; : 


the University of Chicago in 1928, has been active in re- 
search in physical metallurgy and applied x-rays, first at 
the Naval Research Laboratory and now at Carnegie 
Institute of Technology in Pittsburgh. 
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A Mechanical Model for the Motion of Electrons in a Magnetic Field 


ALBERT ROSE 
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(Received June 8, 1940) 


A mechanical model in the form of a gyroscope may be 
made to simulate the path of an electron in a uniform 
magnetic field acted upon by transverse electric fields. In 
using the model, one employs the following substitutions: 
The uniform magnetic field is replaced by the spin velocity 
of the gyroscope, the electric fields are replaced by mag- 
netic fields which give the same field configuration, and 
the electron is replaced by one pole of a permanent magnet 
mounted on the axis of the gyroscope. The spin velocity 
may be arbitrarily chosen for convenient observation. 
Then, depending upon the value of the spin velocity, the 
dimensions of the gyroscope, the strength of the permanent 
magnet and the scale factor used in setting up the model 
the strength of the magnetic field in the model is adjusted 
to correspond to the strength of the electric field in the 
actual case The path described by the magnetic pole, under 


INTRODUCTION 


ECHANICAL analogs of electron motion 

have been found of value not only as an 
aid to pedagogy but also as an aid to original 
investigations in the design of modern electron 
tubes. The use of small balls rolling on a rubber 
membrane to simulate the motion of electrons 
in a two-dimensional electrostatic field is a well- 
known example. While this model has been 
widely used, similar models for electron motion 
in the presence of a magnetic field either have 
not been proposed or have not been convenient 
enough for ready application.! This paper will 
describe a simple model for the motion of 
acted upon by a_ two-dimensional 
electric field in a plane normal to a uniform axial 
magnetic field. The essential part of this model 
and that which traces the electron paths is a 
gyroscope mounted with its axis in the direction 


electrons 


1A mechanical model for electron motion in axially 
symmetric magnetic fields is described by E. Bruche and 
A. Recknagel (Zeits. f. tech. Physik 17, 126-134, 1936). 
This form requires a separate construction of the model for 
each set of initial conditions. 
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these conditions, corresponds to the actual path of the 
electron. Since the error involved in identifying the equa- 
tions of motion of the gyroscope with those of the electron 
is of the order of 6?, where @ is the half-angle of the cone 
through which the spin axis is allowed to move, the range 
of observations must be restricted accordingly. Some 
familiar structures for which electron paths have been 
observed and photographed are deflection plates in a uni- 
form magnetic field, a diode in an axial magnetic field, and 
a split-anode magnetron. By periodically reversing the 
current through the field magnets, it is possible to observe 
the effect of alternating electric fields. In one arrangement, 
in which the field was alternated in resonance with the 
revolution of the electron about the magnetic field lines, 
a path similar to that of charged particles in a cyclotron 
was obtained. 


of the magnetic field.? The first part of the paper 
will compare the equations of motion of an 
electron, in the electric and magnetic fields 
mentioned above, with the equations of motion 
of the end of the axis of the gyroscope acted upon 
by forces normal to its axis. From this compari- 
son, the approximations involved in using the 
model will be set forth as well as the quantitative 
relations between the constants in the two sets 
of equations. The second part of the paper will 
describe the application of the model to some 
familiar types of electron motion and will include 
photographs of the electron paths as traced by 
the gyroscope. 


ANALYSIS OF THE MODEL 


In the cylindrical-coordinate system 7, ¢, 2 of 
Fig. 1, let JZ be a uniform magnetic field in the 


2 It has come to the writer’s attention, since this paper 
was written, that the similarity between the force on an 
electron due to its motion transverse to a magnetic field, 
and the force on a gyroscope due to the angular velocity 
of its axis was suggested in 1933 by Professor E. C. Kemble. 
(See William F. Osgood, Mechanics (Macmillan, 1937), 
p. 234.) 
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Fic. 1. Equations of motion of an electron in a uniform 
magnetic field and of the end point of the axis of a gyro- 
scope. 


z direction, and let E,, Ey, Ez be the electric 
field intensities in the 7, ¢, and z directions. The 
equations of motion of an electron in this 
system are 


*#—r¢?+e/mrdH = —e/m_E,, (1a) * 
rp+2i¢—e/mil = —e/mEs, (1b) 
z= —e/mE.. (1c) 


Of these equations, we shall be primarily 
interested in (la) and (1b) which describe the 
motion of the electron in the r—@ plane, since 
the application of the model is restricted to this 
plane. Fortunately, motion of the electron in 
the z direction is independent both of H and of 
motion in the r and ¢ directions. The z compo- 
nent of the electron path may be obtained from 


* Note: the absolute value of ‘‘e,” the charge on the 
electron, is used throughout this paper. 
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a knowledge of E. and may readily be combined 
with the r—@ plane projection of the path to 
give the complete picture. 

Corresponding to Eqs. (1a) and (1b) for the 
electron, we now desire similar equations for the 
motion of the end point of the gyroscope axis in 
the same r —¢ plane. The motion of the gyroscope 
has been treated extensively both in technical 
journals and in textbooks on mechanics. It will, 
therefore, suffice to indicate here the method for 
obtaining the particular form of the equations 
of motion used in this paper. Lagrange’s equa- 
tions for the general case of a rotating rigid 
body fixed at its center of gravity are 


0 /0T oT 
-(—)- —=]Fy, (2a) 
ot\ 00 06 
0 /o0T oT 
: (- -) ius =L5F;, (2b) 
dt\ dd do 
0/s0T oT 
-( -)-— =LyFy,. (2c) 
Ol\ dy Oy 


The coordinate angles 6, ¢, ¥ are the Eulerian 
angles, two of which, @ and @¢, are shown in 
Fig. 1. The third angle, y, measures the rotation 
of the gyroscope disk with respect to the plane 
formed by the gyroscope axis and the 2 axis 
(Fig. 1). Fe, Fy, Fy, are the forces and Lg, Lz, Ly, 
the moment arms for the three torques. The 
kinetic energy function, J, is most readily 
expressed as 


T= $(Aw,?+Buw.?+ Cw3’), (3) 


where A, B, C are the principal moments of 
inertia and w, we, w3 are the angular velocities 
about the three principal axes of inertia. The 
explicit expression for 07/06, and similar terms 
in Eq. (2), is obtained from 


oT Ow) Owe 0w3 
— = Aw;— + Bux— + Cozs—, etc. (4) 
06 06 06 06 
and from the transformation equations :* 
w,=6sin Y—¢ sin 6 cos y, (5a) 
w2=6cos ¥+¢ sin @ sin y, (5b) 
w3=0+¢ cos 8. (5c) 


3See J. H. Jeans, Theoretical Mechanics (Ginn and Com- 
pany, 1907), p. 347. 
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Fic. 2. Gyroscope set up to record the path of an electron 
in a magnetron. 


These equations, together with the fact that if C 
is the moment of inertia about the gyroscope (or 
spin) axis A is equal to B by symmetry, allow Eq. 
(2) to be written as 
A6+CY¥¢ sin 0 

—(A—C)¢* sin 6cos 0=LeFs, (6a) 


A¢ sin? 6— C6 sin 0 
+2(A —C)¢6 sin 6 cos 0 
+C(¥+¢ cos 0) cos 0=L5F, (6b) 
C(v+¢ cos 6—¢6 sin 6) = LyFy. (6c) 


In the ordinary operation of the gyroscope, if 
one neglects friction, Fy=0 so that 

v+¢ cos 0—¢6 sin 6=0 (7) 
and ae 

¥+¢ cos 6=constant. (8) 


From Eq. (8), the constant is equal to the value 
of ¥ when ¢=0. That is to say, if the gyroscope 
is given an initial spin velocity and then allowed 
to move under the influence of forces transverse 
to the axis, the constant is equal to the initial 
spin velocity. Let W be this spin velocity. Eqs. 
(7) and (8) with the constant replaced by ¥ 
may be used to rewrite Eq. (6) as 

” : LF 

6—¢?* sin 6 cos 0+ C/A¢gW sin 6=—-, (9a) 


4 


“ . a ES : 
¢@ sin 6+2466 cos 0—C/AdV= —" (9b) 


v=0. (9c) 


On the assumption that the point of application 
of the forces Fs and Fy, is at the end of the axis 
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of the gyroscope, Lg and Ly become according 
to Fig. 1, LZ and L sin @. At this point it is well 
to emphasize that no approximations have yet 
been made. Eqs. (9a), (9b), and (9c) are the 
equations of motion of a rigid body, with its 
center of gravity fixed and with its spin axis an 
axis of symmetry, subjected to forces transverse 
to the spin axis—in brief, the equations of motion 
of a gyroscope. If, now, @ is restricted to small 
values such that <1, we can replace cos @ by 
1 and sin @ by r/L where r is the distance the 
end of the gyroscope axis moves, measured along 
the arc of its motion, from its origin at the z axis. 
Since @ is restricted to small values, the ‘‘r’”’ 
defined above may be identified with the r in 
Fig. 1 used as one coordinate of the electron in 
the r—@ plane. With this approximation, (9a) 
and (9b) may be written in a form to be com- 
pared with the electron equations (1a) and (1b), 
namely 





° 
— 


#—rd6?+C/Aro¥ Pn (10a) 


9 


ro+27%¢6—C/ATW <7 Fe (10b) 


Comparison of Eqs. (10) and (1) shows first, 
that the two sets are identical in form and, 
second, that the following correspondences exist 
between their parameters: 


electron gyroscope 

e/mH C/AW, (11a) 
e/mE, L?/AF,, (11b) 
e/mE, L?/A Fy. (11c) 


The formal identity of the two sets of equations 
says, immediately, that if a scale model of the 
force field E,, Ey is set up so that it acts trans- 
verse to the axis of a rotating gyroscope, then 
the path described by the end point of the 
gyroscope is, for small values of @, the path that 
an electron would describe in that force field and 
in a uniform magnetic field normal to the plane 
of the force field. Variations of the spin velocity 
of the gyroscope or of the magnitude or direction 
of parts of the force field will alter the path of 
the gyroscope and allow one to obtain a qualita- 
tive understanding of the changes of path that 
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Fic. 3. Cycloidal path of an electron in crossed electric 
and magnetic fields. 


an electron would undergo for similar variations 
of magnetic and electric fields. This is usually 
the way in which a model is used. It is, never- 
theless, desirable to know also the quantitative 


relations between the “magnetic” and “electric” 
fields of the scale model and the actual magnetic 
and electric fields acting upon the electron. 
A discussion of these relations follows. 

The performance of the model differs from the 
performance of the electron in two significant 
respects. First, the model takes longer to describe 
its path than does the electron. Second, the 
model is usually constructed to operate on a 
larger scale. (Since the model represents only 
one electron, the paths it describes are free from 
the influence of space charge.) The question is: 
“Under what conditions will the electron path 
be unchanged when either (or both) the transit 
time and the scale factor are changed?” By 
“unchanged” is meant either “‘identical’’ when 
only the transit time is changed or “differing 
only by the scale factor’’ when the scale factor 
of the electrodes is changed. These conditions 
may be obtained in the following manner. Let 
the time scale in Eq. (1a) be changed from / to 
t/k where k is an arbitrary constant. Then, 
He/m which has the dimension of 1/7 must be 
multiplied by k and E,e/m which has the 
dimensions of L/7? must be multiplied by & in 
order to leave the solution for the path un- 
changed. In the same way, the conditions may 
be found for leaving the path unchanged when 
the electrode dimensions are altered. The results 
of this analysis may be summed up by the 
statement that the electron path will remain 
unchanged under any of the following conditions. 
(1) Changes of E and H which do not change 
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E/H?. (2) Changes of scale factor which include 
multiplication of E by the same factor. (3) 
Changes of scale factor which include division 
of H by the square root of the scale factor. 
It is to be understood in the above conditions 
that whenever the electrode potentials are 
changed the initial energy of the electron is to 
be changed by the same factor. 

The procedure for setting up a model to 
represent quantitatively a given electron prob- 





Fic. 4. Same as Fig. 3 with the addition 
of an initial velocity. 


lem may now be briefly stated. First, let the 
model have the same linear dimensions as the 
original. Also, let H, E, A, C, ¥, and L be known 
quantities. There remains the adjustment of the 
values of F, and F, so that the gyroscope will 
describe the same path as the electron. To do 
this, the difference in time scale between the 
model and the original is determined by the 
relation (see also (11a)) 


e/mIl=kC/AW, (12) 


which defines k. From (12) and the first condition 
of invariance we can write (see also Eqs. (11a) 
and (11c)) 

e/mE,=k*L?/A F,, (13a) 


e/mE,=k*L?/A Fs. (13b) 


These equations define F, and Fs and permit 
setting up a quantitative model having the same 
dimensions as the original. If the model is to be 
enlarged by a scale factor S over the original, 
then, according to the second invariance condi- 
tion, the forces F, and F, are to be multiplied 
by S. To summarize: 

A model S times as large as the original may 
be set up so that the gyroscope describes corre- 
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spondingly the same path as the electron if limits within which the model may be used, and 
obtained the quantitative relations between the 


PF, = 3 a model and the original. This much would be of 

k? Am academic interest only, unless some convenient 

and or means of setting up a force field to act on the 
F,= . “ne - gyroscope were available. The force field must, 

k? Am of course, approximate the electrostatic field in 

where the original. Fortunately, such a field can easily 
ben ” - be set up for many types of electrostatic fields 


Cm ¥ by means of electromagnets with specially shaped 
pole faces. For example, to represent the electro- 
APPLICATION OF THE MODEL 


The first part of this paper has described a 
model for tracing electron paths, defined the 





Fic. 7. Path of an electron in a magnetron. The lower plate 
; be I 
is the less positive plate. 





Fic. 5. Path of an electron in a central force field. 





Fic. 8. Path of an electron in a magnetron showing 
electron ‘‘trapped” for an appreciable time in a low electric 
Fic. 6. Same as Fig. 5 with large initial velocity. field region between cathode and the lower potential plate. 
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Fic. 9. Path of an electron beam in a magnetron as 
obtained by G. R. Kilgore using gas to make the beam 
visible. 


static field in a diode, an electromagnet with 
one pole connected to a soft-iron cylinder and 
the other to a soft-iron rod along the axis of the 
cylinder has been used. The electron has been 
represented by one pole of a long thin permanent 
magnet mounted so that one end is at the end 
of the gyroscope axis and the other end near 
the center of the gyroscope. This end contributes 
comparatively little to the moments of force 
acting on the axis both because it has a short 
lever arm and because it is somewhat removed 
from the external magnetic fields. 

Another problem of practical importance is to 
mount the gyroscope so that there is a minimum 
of friction in the bearings. Friction in the spin 
axis bearings causes a slowing down of the spin 
velocity or in terms of the original, a continu- 
ously decreasing magnetic field. Friction in the 
gimbal bearings tends to damp out any rapid 
oscillations of the gyroscope and, in larger 
amounts, tends to introduce objectionable ex- 
traneous forces which permit the gyroscope axis 
to drift in the direction of the applied forces. 

The gyroscope used to obtain the following 
paths is the small model made by the Sperry 
Gyroscope Company. The spin axis bearings 
were found to have sufficiently little friction so 
that even for small spin velocities the ‘‘magnetic 
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Fic. 10. Path of an electron in crossed electric and mag- 
netic fields when the electric field is alternating in resonance 
with the angular period of the electron. 


field’’ remained substantially constant through- 
out a run. It was found necessary, however, to 
provide ball-bearing mountings for the gimbals 
to reduce their friction to a usable amount. The 
and field apparatus to 
represent a are shown in 


gyroscope associated 
magnetron Fig. 2. 
The end of the axis toward the ‘‘magnetron”’ 
has mounted on it a long thin Alnico magnet, 
the outer end of which represents the electron. 
A small flashlight bulb and battery are mounted 
on the other end of the axis for the purpose of 
photographing the path. A light chopper in 
front of the camera provides the time axis in the 
recorded picture of the path. 

Figure 3 shows the familiar cycloidal motion 
of an electron in crossed electric and magnetic 
fields. In Fig. 4, the same type of motion is 
recorded except that the “‘electron’’ was given 
an initial velocity. In all of these photographs 
the uniform magnetic field is normal to the plane 
of the picture. 

The path of an electron, started from rest, in 
the radial electric field of a diode is shown in 
Fig. 5. In Fig. 6, the electric field is radial but 
generated by a point charge rather than a line 
charge. The electron was started very near the 
point charge and given an appreciable initial 





4This type of motion is to be found in the electron 
multiplier described by V. K. Zworykin, G. A. Morton, 
and L. Malter, Proc. I. R. E. 24, 361-370 (1936), and in 
the mass spectrograph described by W. Bleakney and 
J. A. Hipple, Jr., Phys. Rev. 53, 521-529 (1938). 
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velocity. The gradual reduction in amplitude of 
successive loops is due to the small residual 
friction in the gimbal bearings. Two typical 
electron paths, under static conditions in a 
magnetron are shown in Figs. 7 and 8. The 
arrangement which generated the magnetron 
field can be seen in Fig. 2. The two plate sections 
were connected to opposite poles of an electro- 
magnet, the upper plate being the more positive. 
A separate electromagnet in the form of a soft- 
iron rod served as the cathode. In both figures, 
it is to be noted that the electron lands on the 
lower potential plate. Figure 8, in particular, 
shows the electron trapped for an appreciable 
time in a region of low electric field between the 
cathode and lower potential plate. For com- 
parison, an actual electron path in a magnetron 
as traced by an electron beam in gas is given in 
Fig. 9.5 Figure 10 is of particular interest in 
showing the ability of the model to reproduce 
electron paths in alternating electric fields. A 
uniform magnetic field supplied by an electro- 


5 This figure is reprinted by courtesy of the Institute of 
Radio Engineers from the paper, ‘‘Magnetron oscillators 
for the generation of frequencies between 300 and 600 
megacycles,”’ by G. R. Kilgore, Proc. 1. R. E.24, 1145 (1936). 





magnet was periodically reversed by means of a 
motor-driven commutator in resonance with the 
motion of the gyroscope. Starting from rest, 
then, the electron described successively large 
circles as it gained energy from the alternating 
field. The fact that the light streaks are approxi- 
mately lined up radially indicates that the 
angular velocity remained constant as the linear 
velocity was increased. (The angular velocity of 
an electron describing circles in a magnetic 
field is determined only by the magnetic field 
strength.) The slight ellipticity of the paths is the 
result of not having the moments of inertia of 
the gyroscope about the two sets of gimbal 
bearings exactly the same. The slight drift of 
the center of revolution towards one side is the 
result of a small d.c. electric field superposed on 
the alternating electric field. 

These few illustrations show that the gyroscope 
may readily be used to give a semiquantitative 
demonstration of various types of electron paths 
in a uniform magnetic field and that for each 
type the initial velocities and field strengths 
may be conveniently varied to explore a wide 
range of conditions including zero magnetic field 
and alternating electric fields. 





Notes on the Power Factor Testing of Long Lengths of Cables 


Eric A. WALKER 
Department of Electrical Engineering, Tufts College Engineering School, Tufts College, Massachusetts 
(Received February 16, 1940) 


The ratio of watts to volt-amperes does not give the power factor of long cables because of 
the resistance of the conductor and sheath. The error involved varies as the length squared 
and might be quite appreciable at lengths of 10,000 feet. Correction factors are given involving 
the use of an infinite series for exact results and a simple algebraic formula for approximate 


results. 


INTRODUCTION 


NUMBER of tests have been devised for 
determining the condition of insulation. 
The best known of these tests are the following: 
first, a direct-current resistance test made by 
means of a megger at several hundred volts;' 
second, a power-factor test made by means of a 
bridge or a watt-meter set ; lastly, an overvoltage 
1 Megger Insulation-testing Instruments (J. G. Biddle 
Company, Philadelphia, Pa., 1928). 


e 
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test. In this test a voltage several times the 
rated voltage of the equipment is applied to the 
insulation for a predetermined period of time. 
During this test the resistance may also be 
measured. 


THE ADVANTAGES AND DISADVANTAGES OF THE 
POWER FACTOR TEST 


There are a number of reasons why the 
power-factor test is preferred. First of all, this 
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Fic. 1. Apparent power factor of a theoretical no-loss 
1,000,000-cm 5-kv cable at 60 c.p.s. 


type of test is made to determine the breakdown 
strength of the insulation and to ascertain its 
probable useful life. There appears to be little 
or no correlation between either of these prop- 
erties and the direct-current resistance. However, 
both the life and the breakdown strength have 
been correlated with the power-factor of some 
materials.2~* Also power-factor is known to show 
certain types of faults which direct-current 
resistance tests will not reveal. The power-factor 
is independent of the amount of material being 
tested. As the power-factor test can be conducted 
at any voltage which the capacity of the testing 
equipment will allow, a voltage which will not 
injure nor destroy the insulation can be chosen. 
Lastly, the power-factor test has been used 
successfully on bushings and transformers in the 
field,» and on short lengths of cable in the 
laboratory.®7? A logical extension of the method 
is to apply the test to longer lengths of installed 
cable and to cables which have been in service 
in order to determine their present condition and 
life expectancy. 

The limitations of the power-factor test are 
obvious. It is a test which gives a measure of 
the average condition of the cable. In other 
words, the test will not reveal a small section of 





?C. L. Dawes and F. L. Hoover, ‘‘Ionization studies in 
paper insulated cables. I,’’ Trans. A.I.E.E. 45, 141 (1926). 

§C. L. Dawes, H. H. Reichard and P. H. Humphries, 
“Ionization studies in paper insulated cables. II,’’ Trans. 
A.1.E.E. 48, 382 (1929). 

4B. L. Goodlet, F. S. Edwards and F. R. Perry, ‘‘Dielec- 
tric phenomena at high voltages,” J. Inst. Elec. Eng. 69, 
695 (1931). 

5]. W. Gross, “Field testing of bushing and transformer 
insulation by the power factor method,” Trans. A.1.E.E. 
57, 589 (1938). 

6 G. B. Shanklin, ‘‘Compensated dynamometer wattmeter 
method of measuring dielectric energy loss,’’ Gen. Elec. 
Rev. 19, 842 (1916). 

7E. S. Lee, ‘The use of a dynamometer wattmeter,” 
Trans. A.I.E.E. 45, 746 (1926). 
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deteriorated cable if it is tested as a part of a 
longer cable which is almost wholly in a satis- 
factory condition. In this respect the test is 
inferior to the overvoltage test which gives the 
breakdown strength of the weakest point if 
failure does take place. The second limitation is 
in the combination of high kva capacity which 
must be provided and the small amount of actual 
power which must be measured. For instance, 
a 5-kv varnished cambric cable with 1,000,000- 
em conductor and insulation 0.109 inch thick 
will require about 4.25 kva charging capacity 
per 1000 feet of length if the test is made at 
5 kv and 60 c.p.s. while the watts loss in the 
insulation at 0.02 power-factor will be of the 
order of 85 watts. However, apparatus does 
exist and tests have been made on lengths of 
cable up to five miles. 

The third limitation, and the one with which 
this paper is concerned, is the effect of the 
resistance and inductance in the circuit being 
measured. 


THE Circuit CONSTANTS 


In testing a cable by power-factor methods the 
path taken by the charging current must be 
studied in order to evaluate the circuit constants. 
Usually the charging current is fed from one 
end of the cable. This current flows out along 
the conductor, and returns through a path 
determined by several factors. If the cable is 
isolated and lies in a fiber duct the return path 
will be through the lead sheath. This makes the 
resistance of the return path quite appreciable. 
Of course, the resistance of the total return 
path may be much lower than the resistance of 
the lead sheath because of neighboring cables or 
other structures which provide an alternative 
ground return, or because of a copper binding 


Apparent Power F actor 





Length in Thousands of Feet 


Fic. 2. Apparent power factor of a theoretical no-loss 
0000 50-kv cable at 60 c.p.s. 
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tape around the insulation of the cable. One can 
only say then that the total resistance cannot be 
lower than the resistance of the center conductor 
or higher than the resistance of the conductor 
and the lead sheath. This resistance can be given 
in ohms per unit length or in ohms for the entire 
structure. For instance, the resistance per 1000 
feet of conductor of the above cable is about 
0.387 ohm. 

The inductance of the circuit, if the current 
returns through the sheath, can be computed as 
the inductance of a concentric cable. An approxi- 
mate formula, which neglects the thickness of 
the sheath is: 


R, 
1= 0.0305] 0.5-+2 In | millihenries per 
R, 1000 feet. 


If the return is not through the lead sheath but 
through other cables, or through ground, the 
inductance will be higher. The above formula 
gives the lower limit only. It is fortunate that 
later developments will show that this factor is 
not important. 

The conductance and the susceptance of the 
insulation of the cable are the characteristics of 
the cable it is desired to measure. For small 
defect angles the power-factor is given by the 
ratio of these quantities: 


: Conductance G 
Power-factor =-—---—-—- = — (1) 
Susceptance B 


The question is, ‘‘Does a measurement at the 
end of a cable give this ratio?” In an attempt 
to answer this question the cable will be con- 
sidered as an evenly distributed impedance, 
commonly called a smooth line. 


THE DEVELOPMENT OF AN EXACT AND AN 
APPROXIMATE FORMULA 


Any measurement at the end of a cable will 
measure the input impedance of the cable, 
whether it is determined by a bridge or by a 
measurement of volts, amperes, and watts. The 
input impedance of a smooth line with the far 
end open is 





Zo 
Siegut — ’ (2) 
tanh 6 
or 
Y input = Yo tanh 8, (3) 
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where 


Zo=(Z/Y)*; Yo=(¥/Z)!; 0=(ZY)}, 


where Z=R+jLw=R+jX, the total series im- 
pedance of the cable, and Y=G+jCw=G+JB, 
the total parallel admittance of the cable. 

Equation (3) can be expressed in a more 
convenient form: 


tanh 6 
Vinput= Yo tanh 6= Y——. (4) 
6 


This is an exact equation for any smooth line. 
The equation can be expanded in an infinite 
series : 


ZY 2 
Y input = v1 ———_+—7Z°?Y?. . | 
- 
(RG—BX)+j(RB+GX) © 
~ccxnf omens 





A test with actual numerical values shows that 
for a wide range of values this equation reduces 
to the following approximate formula: 


RB? 
Y input = G+— ~ -+7B. (6) 
Then Eq. (1) becomes: 
G RB 
Observed power-factor =—-+—. (7) 
3 


Let the observed power factor = 6’ and the real 
power factor G/B=6. Then 


RB 
§=5’-——_. (8) 
3 


This equation shows that the observed power- 
factor is too high by a correction factor depending 
on the resistance and susceptance of the cable. 
It does not depend on the inductance. The 
correction factor also increases as the square of 
the length of the cable. 

If the volts amperes and watts are read for 
the cable the equation becomes: 


9 


watts I?R 
6 =—— {1- | (9) 


volt-amperes 





In Fig. 1 the apparent power-factor of the cable 
already described is shown for lengths up to 
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100,000 feet, on the assumption that the cable 
had zero power-factor. This figure then shows 
the amount which must be subtracted from the 
observed power-factor of a cable of these 
dimensions in order to determine the real 
power-factor. In this figure, the correction factor 
as determined by Eq. (8) is indicated by the 
crosses, while the circles give the exact correction 
factor found by using as many terms as are 
significant in Eq. (5). It is seen that the agree- 
ment is fairly good. 

Figure 2 shows the same sort of curve for 
another type of cable. This cable has dimensions 
quite different from those of Fig. 1 and it was 





CONCLUSIONS 


For long cables the power-factor as observed 
at one end must be corrected to give the true 
power-factor of the insulation. 

Although the correction factor is not important 
for short lengths, it is important for lengths of 
cable of the order of 10,000 ft. 

The resistance of the return path of the current 
has considerable influence on the correction 
factor. 

The inductance of the cable has little effect on 
the measured values. 

The correction factor will be higher at higher 
frequencies. 


chosen to show that Eq. (8) is valid for widely 


A simplified equation for computing this 
different types of cables. 


correction factor is quite satisfactory. 





Similarity of the Stress Distributions in a Circular Disk and a Square Plate 


H. Oxuso 
Technical College, Yonezawa, Japan 
(Received April 18, 1940) 


_—. a solution has been carried out for a circular disk compressed diametrically. But few 
solutions for a square plate were found formerly. In this paper we will introduce a solution 
which is more convenient for numerical calculations than the solutions of the previous authors, and 
show that the stress distribution in a square plate compressed by a pair of opposite forces acting on its 
two sides is similar to the one in a circular disk compressed diametrically. 

The generalized plane stress in a plate is generally given by! 


X.=—lim @[k*ei(x+iky) + ¢2(x+iy)], 


Y,=lim RL gi(x+iky) + go(x+iy) ], (1) 
X,=lim Tl Roi(x+iky) + ¢o(x+iy) |, 


where ¢g; and ¢2 are two arbitrary functions determined to satisfy the boundary conditions, ® and T 


denote, respectively, the real part and the imaginary part divided by 7 of a function directly following 
them. 


When a square plate of side length 2 is compressed by a pair of opposite forces P acting on its two 
sides at their centers as is shown in Fig. 1, the boundary conditions are 


Xs 


0, X,=0, when x=+1, 


P “ (2) 
~~ P 2 cw ane, X,=0, when y=+1. 


Y, 


n=1 


1H. Okubo: Science Reports of the Téhoku Imperial University, Series 1, XXV. No. 5 (March 1937). 
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Stresses on the axes 


Cc 
» 





TABLE I. 
Po = —0.5 
P,; =—1.047 QQ: = 0.255 
P, =—0.963 OQ» =—0.116 
P; =—1.026 QO; = 0.051 
P, =—0.981 O, = —0.029 
P, = —0.989 O. =—0.010 
P; =—1.009 Q; = 0.005 
P, =—0.993 Os = 0.000 
P, =—1.006 Oy =—0.002 
Po= —0.996 Oro= —0.001 
Pi,;=—1.004 Q11= —0.006 
P..= —0 997 VO2= 0.007 
= = ———_— _ Fic. 2. 
Now putting 
0 7 
gi(x+iky) => A, cos nx(x+iky), 
n=0 
wo (3) 
go(x+iy) => B, cos nr(x+7y). 
n=C 
The stresses (1) satisfy the boundary conditions 
@ 
Y,=)> Pacosnrx, X,=0, when y=+1, (4) 
if n=0 
P,, sinh nx ) 
cosh urk sinh nx—k sinh nk cosh nr (5) 
> « 


_ —kP,, sinh nrk 
*” cosh nnk sinh nr—k sinh nk cosh nn J 
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TABLE II. 





x y Xx Vy 

0 1.00 — 2 

0 0.75 0.328 — 2.551 
0 0.50 0.298 — 1.322 
0 0.25 0.296 — 1.017 
0 0 0.296 — (0.939 
0.25 0 0.231 — 0.800 
0.50 0 0.103 — 0.503 
0.75 0 0.027 —(.214 
1.00 0 0 0 


(5) into (1), we have 


® 
X,=)>,. P,&,(y) cos nrx, 


n=0 


Inserting A,, B, given by 


=>> P.in,(y) cos nx, > 
n=O 


4) 
=> P, f(y) sin nxx, 
n=0 

where 


E,(y) =— —_ ~te 


sinh nz porn aehae 


xe vy" 
0. 318 — 2 
0.318 — 2.592 
0.318 | — 1.379 
0.318 — 1.040 
0.318 —0.955 
0.248 — 0.809 
0.115 — 0.497 
0.025 — 0.202 
0 0 


sinh nz cosh nxyt+nny sinh mz sinh nry—nz cosh nr cosh nry 


sinh mz cosh nry— ry sinh nx sinh nry+nr cosh nm cosh nry 


"n(y) =————- — ——— 
rary nT chal oetas 


nx cosh nr sinh nry—nry sinh nr cosh nny 
¢.(y) =———_ 





sinh mz cosh uxr+nur 





On the sides x = +1, the shear vanishes and the normal tension is 


a z= P,(- 1)"é,(y), 


n=l1 


which is expanded in Fourier’s cosine series and given by 





=-E Q;cos jny, 
j=1 
where ee 
~% 4(—1 1)‘nj? sinh? nx 
=o, Pal—1)ridte, Dp, I, j)=——— —_—___—__—_—— 
a(n*+ 77)? )2(sinh nx cosh ax-+nx) 
Similarly we can find the stresses which satisfy the following boundary conditions 
X,=> QO; cosjry, X,=0, whenx=+1, 
, as j=l 
= x Qjn;(x) cos jry, 
j=l 
Y,= Qj§i(x) cos jry, 
j=1 
yv=D Ojf;(x) sin jry. 
I7=1 J 
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(6) 


(9) 
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On the sides y= +1, the shear vanishes and the normal tension is 


Y,=> 0)(—1)§(x) = ¥ Rn cos mrx, 
j=1 m=1 


where 
Rn => Qi(—1)47(Gj, m). 
7=1 
The stress system (6) superposed with (9) satisfies the boundary condition (2), if 


P,t+ > P,,(—1)"*} > (—1)‘J(Gj, n)I(m, 7) +P =0. (10) 
j=1 


m=lL 


From Eq. (10) we can determine P,,, and inserting them into (8) we can find Qj, all the unknown 
constants are known. Values of P, and Q; calculated from Eqs. (8) and (10) are given in Table I. 
The stress system (6) superposed with (9) is 


«x 


X;= | PrEn(y) cos ntx+Qyn»(x) Cos ny}, | 
n=1 
P xr 
Y,=——+) {Pruin,(y) cos ntx+Q,£&,(x) cos nry}, ! (11) 
Z n=1 
x 
X,= > {Pafn(y) sin nxrx+Q,¢,(x) sin ny}. | 
n=1 J 


Values of stresses on the x and y axes are calculated from the above equations and tabulated in 
Table II and are plotted in Fig. 2 (* shows the values in a circular disk). From the figure it is obvious 
that the stress distribution on the axes are closely similar with the one on the horizontal and vertical 
diameters of a circular disk of radius unity under the cognate boundary conditions. If that is com- 
pared with the fact that the stress distribution in an elliptic plate with an eccentricity near unity is 
not similar to the one in a circular disk,” the result is very interesting. 

In conclusion, the writer wishes to express his thanks to Mr. S. Oba, Head of the Yonezawa 
Technical College for his encouragement throughout the investigation. 


2H. Okubo, Phil. Mag. 27, 508 (1939). 
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Stress Analysis by Three-Dimensional Photoelastic Methods* 


DaniEL C. DRUCKER AND RAymMonp D. MINDLIN 
Department of Civil Engineering, Columbia University, New York, New York 
(Received June 13, 1940) 


This paper describes the results of an investigation and 
extension of the various proposed procedures and methods 
of analysis for the photoelastic determination of three- 
dimensional states of stress. Two important limitations of 
the previous developments are removed. A true extension 
to three dimensions is given by the determination of the 
effect on retardation of appreciable variation in the 
orientation of the secondary principal stresses in planes 
perpendicular to the wave normal (such variation being 


INTRODUCTION 


WO general practical methods have been 
proposed, and to some extent used, for the 
photoelastic determination of three-dimensional 
states of stress. In one, the stresses are “‘frozen”’ 
in the model and slices cut from it are analyzed; 
in the other, the phenomenon of scattering of 
light is employed. The technique of measure- 
ment required is different for the two cases and 
therefore they will be discussed separately. How- 
ever, both methods depend upon the effects of 
stress on the propagation of light, and these are 
the same regardless of the method of observation. 
In the previous developments of the theory, only 
those effects were considered which are present in 
ordinary two-dimensional states of stress. 

These effects in three-dimensional states of 
stress differ considerably from those present in 
the usual two-dimensional case. The rotation 
along the wave normal of the axes of secondary 
principal stress' in planes perpendicular to the 
wave normal.is the principal cause of the differ- 
ence and has not been considered previously. 

Knowledge of what occurs in general states of 
stress permits broadening of the methods of ob- 
servation to include entire planes rather than 
confining attention to a single point at a time. 


WELLER’s METHOD 


If one wishes to analyze a three-dimensional 
state of stress in a model under actual load, the 


* Publication assisted by the Ernest Kempton Adams 


Fund for Physical Research of Columbia University. 
1R. D. Mindlin, ‘A review of the photoelastic method of 
stress analysis,” J. App. Phys. 10, 229 (1939). 
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the rule rather than the exception for general states of 
stress). Also, methods for the analysis of whole planes are 
presented, thus avoiding the cumbersome and tedious 
point by point procedures that have been advanced. In 
the particular case of plane stress, these methods reduce 
to a purely optical technique for determining the principal 
stresses themselves, or more easily their sum, for the entire 


model. 


only practical method of general application 
available at present is that proposed by Weller.? 
This method is based upon the fact that in a 
transparent medium, light scattered at right 
angles to the incident wave normal is plane 
polarized, the amplitude of the vibration reach- 
ing the eye of an observer being proportional to 
the projection of the light vector on the plane 
perpendicular to the direction of observation. 
The result is similar to that produced by the 
analyzer in a polariscope. (See Fig. 1.) If the 
light incident at a point is plane polarized, the 
intensity of scattered light observed at right 
angles to the wave normal will thus vary from 
zero when looking along the direction of vibration 
to a maximum when looking perpendicularly to 
this direction. 

When photoelastic materials are stressed they 
become doubly refracting and, for a given wave 
normal, light can be transmitted only in vibra- 
tions along two mutually perpendicular axes 
lying in the planes perpendicular to the wave 
normal. These directions, which we may desig- 
nate as the axes of secondary principal stress! 
are determined in the same manner as the 
principal directions for plane stress. (The normal 
and shearing stresses acting on the planes per- 
pendicular to the wave normal need not be 
considered in measuring or calculating secondary 
principal stresses as the former may be varied 


2 R. Weller, ‘‘A new method for photoelasticity in three 
dimensions,” J. App. Phys. 10, 266 (1939). ‘‘Three- 
dimensional photoelastic analysis by scattered light,” 
Papers presented at the 9th meeting of the Eastern Photo- 
elasticity Conference, May 1939, p. 19. 
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at will without changing the latter, even though 
such variation changes the magnitudes and direc- 
tions of the true principal stresses.) The light 
scattered from any point (O, Fig. 1) can there- 
fore be considered as coming from two plane 
polarized vibrations (along p and q) differing in 
phase by an amount depending upon the in- 
tegrated optical properties of the material along 
the path from the point of entrance (O’) into the 
model to the point at which the scattering 
occurs. In general the two vibrations combine to 
form an elliptical vibration, the plane of the 
ellipse being perpendicular to the wave normal. 
The intensity of the light scattered at right 
angles to the wave normal is proportional to the 
square of the length of the projection (e.g., QQ 
or 77) of this ellipse on a plane perpendicular to 
the direction of observation. 

If a whole plane of the model is illuminated by 
a sheet of parallel rays of polarized light formed 
by passage through a polarizer, collimator, and 
a slit, the light scattered in a given direction 
perpendicular to the wave normal will in general 
vary in intensity from point to point in the 
plane depending upon the projections of the 
ellipse of vibration associated with each point in 
the plane, the alternate lighter and darker 
regions forming a fringe pattern. This inter- 
ference pattern will vary in distinctness and even 
in form as the direction of observation is changed, 
but in general it may be said that the distance 
between fringes in the direction of travel of the 
light is a function of the state of stress between 
the points of minimum intensity. 

Weller’s method for a general state of stress is 
to fix attention on any point desired and rotate 
the model until the minimum distance between 
fringes straddling the point is obtained. If the 
stress field is perfectly uniform in the vicinity, 
this minimum distance gives the maximum prin- 
cipal stress difference at the point according to 
the simple relation’ 


g3—~01>~-, (1) 


*R. Weller and J. K. Bussey, ‘‘Photoelastic analysis of 


three-dimensional stress systems using scattered light,” 
lech. Note No. 737 N. A.C. A. 


VOLUME 11, NOVEMBER, 1940 





where o3 and go, are the algebraically largest and 
least principal stresses, respectively, C is the 
stress optical coefficient in stress units per unit 
thickness per fringe, and d is the distance be- 
tween fringes measured along the wave normal, 
i.e., along the direction of the intermediate 
principal stress o2. The formula and reasoning 
are exactly the same as in two-dimensional work. 

If os and o; vary in magnitude along the wave 
normal, but remain fixed in direction, the stress 























Fic. 1. Scattering of light. The entering vibration OJ’ 
breaks up into two components, one along p, the other 
along g. At any point O in the model the two components 
in general will combine to form an elliptic vibration with 
axes Oi and Oj in the plane perpendicular to the wave 
normal O’O. The intensity of scattered light seen by an 
observer is proportional to the square of the projection of 
the ellipse (e.g., QQ or TT) on a plane perpendicular to 
the direction of observation. 


found from Eq. (1) will exist at some point inter- 
mediate between the two fringes. For small 
gradients of stress magnitude, the midpoint is 
close enough. If greater accuracy should be 
desired, plotting the fringe positions against their 
order will determine the point more closely. 

If the orientations of the principal stresses do 
not vary, there is no theoretical objection to 
Weller’s proposal. However, as will be shown in 
the following section, if there is rotation along 
the wave normal of the axes of primary or 
secondary principal stress in planes perpendicular 
to the wave normal, the distance between fringes 
is no longer dependent solely on the magnitudes 
of the stressés and Eq. (1) or its equivalent for 


secondary principal stress is no longer. valid. 
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EFFECT OF ROTATION OF AXES ON RETARDATION 


As we are interested in determining the effect 
of rotation of the axes of primary or secondary 
principal stress, it is desirable to strip the 
problem of all stress variables except orientation. 
The fundamental case to be considered is thus 
one in which the magnitudes of the stresses re- 
main fixed while their orientations, in planes 
perpendicular to the wave normal, change at a 
constant rate along the wave normal. In any 
case we have to solve a problem of wave propaga- 
tion in a nonhomogeneous, nonisotropic medium. 
The mathematical formulation of such problems 
generally leads to a system of simultaneous, 
linear, partial differential equations of the second 
order with variable coefficients. 

The differential equations may be formulated 
in terms of an elastic ether theory of light. 
From a general point of view, this theory is not 
the most satisfactory one available; but it is 
very convenient for photoelasticity as it is based 
upon the ideas of stress, strain, and displacement, 
and upon the equations of motion with which 
workers in elasticity and photoelasticity are 
familiar. The resulting differential equations are 
exactly the same as those obtained from the 
electromagnetic theory of light if the symbols 
representing displacement, density, and shear 
modulus are interpreted as magnetic force, mag- 
netic permeability, and reciprocal of dielectric 
coefficient, respectively. 

The equations of motion in terms of the stresses 
are 


Oo,/ OT, Oz,’ 07u' 

ei aes 

Ox Ov 0s ot |} 
Otry Oa, OTyz' 0°" 

=p ‘a (2) 

Ox ov Os or 
Orz:;, OTy: 9a,’ 07m’ | 

ne ar 

Ox ov Os ot? J 


where a,’, oy’, o2', Try’, Tyz, Tr: are the ether 
stresses, u’, v’, w’ are the ether displacements, 
and p is the density. 

The ether stresses may be expressed in terms 
of the ether displacements and the elastic 
coefficients of the ether by three equations of the 
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type’ 


dw’ du’ du’ = av’ 
—E; + —E,{ —-+——}. (4) 
Ox dz Ov Ox 
The elastic coefficients (the E; and G;) are 
variables for a nonhomogeneous medium. They 
are related to the principal wave velocities a, b, c 
by equations of the type 
Gy’ 
a°= ‘ (9) 
p 


in which Gy’ is a stationary value of the G;. 

Finally, the principal velocities a, 6, c are 
determined from the stress-optical laws’ by 
three equations of the type 


a? =v9" +a(o2+03) +o, (6) 


in which vp is the velocity in the unstressed 
medium and @ and 8 are stress-optical coefh- 
cients, constant for a given material and wave- 
length, A, of light. 

Considering a Cartesian coordinate system 
x, y, with the origin at entrance to the medium, 
the simplest state of stress that will satisfy the 
conditions of the present problem is a simple 
tension, ¢=03, in the yz-plane inclined to the 
z axis at a variable angle g=(A/2)x, where A is 
a constant. If the wave normal is in the direction 
of the x axis, A=(21r0/C)x=2kx is the relative 
retardation in radians caused by the magnitude 
of o alone. The ratio of rotation to retardation 
is ¢/A aconstant which we shall call R/ 2. 

Using these specifications for the state of 
stress, and combining Eqs. (2), (3), (4), (5), and 
(6) we find the following simultaneous differential 
equations for v’ and w’: 


4M. G. Lamé, Legons sur la Théorie Mathématique de 
l’Elasticité (Paris, 1866), p. 231. 

5 E.G. Coker and L. N.G. Filon, Treatise on Photoelasticity 
(Cambridge University Press, 1931), p. 198. 
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02y’ 0°y’ 0°y’ Ov’ 
= /]—_— B cos Ax ——+ BA sin Ax 
or Ox? Ox? Ox 
0?w’ Ow’ 
+B sin Ax——+BA cos Ax— 
Ox? Ox 
(7) 
07w’ 02w’ 0°w" Ow’ 
— = /]—_- + B cos Ax —_ — BA sin Ax—— 
ot? Ox? Ox? Ox 
Oy’ dv’ 
+B sin Ax——-+BA cos Ax— 
Ox? Ox) 
in which 


a—p a+p 
B( _— Je I =v+( reas Je 
2 2 


Equations (7) may be reduced to equations 
with constant coefficients by referring the de- 
pendent variables to the axes of principal stress. 
Calling the displacements along the 2 and 3 axes 
(perpendicular and parallel to the direction of o) 
v and w, we find 


d*vy A?D Ow dz ) 
F = v+HA— =— 
Ox? 4 Ox = or? 
(8) 
0*>w <A?*F dv Ow | 
felis taal | 
Ox? 4 Ox ot? 


in which 


F=0,)°+60=H—B; D=v,)?+ac=11+B. 


The independent variables may be separated 
in the usual manner by setting 


and, finally, the dependent variables V and W 
may be separated by differentiation and addition 
of the two resulting equations. This results in 
the same fourth-order linear differential equation 
with constant coefficients for both V and W: 
atV A? a7V 
D F—-+]| (D+ F)w? -—(D?+ F*) + IPA? |— 
ox' 4 Ox? 


v= V(x)et! 


w= W(x)eri+t 


A*w? 44 


+0 —~—(D+F)+- pF =0. (10) 
4 16 


The exact solution of this equation is so 
cumbersome that it is quite unusable. It can 
however be simplified by taking into account the 
orders of magnitude of the quantities involved. 
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If the product 2k is much less than unity, most 
of the terms in the solution disappear. Even in 
highly stressed Bakelite, 2kX will be less than 0.01 
so that the simplification is justified, especially 
so as the square of the term appears and not the 
term itself. 

Considering only the forward traveling waves 
and taking v= M cos wt at x=0, we then find 


Ww 
v=[K sin kSx+M cos kSx | cos («- ~«) 
VH 


+(L sin kSx) sin («- <x), (11) 
VH 
where S=(1+R?)? and K and L are constants. 
Taking w= WN cos wt at x=0, the same equation 
will hold for w except for constants. Thus there 
are four constants to be determined. They can 
be obtained in terms of M and WN by substi- 
tuting in Eqs. (8). When this is done, apparently 
eight relations exist among the four constants. 
Making the assumption above as to the magni- 
tude of the retardation and also a similar one 
for the magnitude of A=2¢/x and also for 
S=(1+R?)!, the eight relations reduce to four 
which determine the constants uniquely. Even if 
a rotation of axes of 360° takes place within a 
distance as small as a millimeter, and if 2k is as 
small as one-tenth of a fringe per inch, \A and AS 
are small compared to unity and can be neglected, 
so that the assumptions are entirely justified. 
v and w can then be expressed in terms of their 
initial amplitudes M and N as follows: 


R 
o=| ——WN sin kSx+ M cos bss 
S : 


w 
xX cos ( _ —x) 
/H 


M w 
_ ( — sin ksx) sin («-—x) 
> JIT 
(12) 





R 
w= +—M sin kSx+WN cos ks] 
S 


w 
Xcos («- = —x) 
VH 
N w 
+(= sin kSr) sin («-——) 
S VHF} 
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assuming them to be in phase on entrance. If ¢ 
is set equal to zero these reduce to the expressions 
which would be obtained for simple tension 
along z. If v and w have an arbitrary phase 
relation on entrance, it is only necessary to add 
a phase constant ¢ to the terms multiplied by / 
and another phase constant e’ to the terms 
multiplied by N. For example, if circularly 
polarized light enters the model, 7,2 is added to 
the phase of either the M or N terms depending 
upon whether the polarization is right or left, 
and M is then set equal to NV. The term (w/\/ H)x 
appears in all of the phase angles with the same 
sign and therefore has no the 
relative retardation. 

From Eqs. (12) for the components of vibra- 
tion, or more clearly from the combined form 


R a ) 
= ——WN sin kSx cos («i : v) 
Ss JH 
: R?2 j 
+] (1 sin RSx) | 
S? | 
w tan kSx 
cos & —x-+arctan ( - —_ )] 
VI S 


R w 
w= +—HWM sin kSx cos («'- - ~«) 
S VH 


1 


R? a 
+[(: —— sin? RSx) 
5S? 
w tan kSx 
X cos |“ —x — arctan ( — )| 
Vil S J 


it may be seen that, wherever tan kSx=0, v and 
w are in phase. This requires kRSx=nz, or 
x=nr/kS. The distance between points in phase 


influence on 





is thus d=2/kS. Now r/k is the spacing one 
would find if only the magnitude of the stresses 
affected the retardation, so that the factor 1/S 
is the effect of rotation. Hence the distance 
between points of equal phase is decreased, or 
conversely the retardation is increased by rota- 
tion. Intensity expressions for either plane or 
circularly polarized light on entrance can be 
set up. From these it is found that minima of 
intensity occur at intervals of r/kS for observa- 
tion in any direction fixed in space or for a 
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direction rotating with the axes of principal 
stress. 

The same results were obtained by a quite 
different method of approach which is of value 
because of its closer physical connection with 
the problem. A series of parallel identical crystal 
plates (corresponding to constant magnitude of 
principal stresses) each with its polarizing axes 
inclined at a constant angle to those of the 
previous plate (uniform rotation of the axes of 
principal stress) was considered. If the number 
of plates is made infinite while the total relative 
retardation due to the magnitude of the stresses 
and the total rotation remain finite, the condi- 
tions of the problem are satisfied. 

There are two essential parts to the solution. 
First, general expressions are obtained for the 
vibrations along each of the two crystal axes in 
any plate by a deductive-inductive process; 
second, the limits of these expressions, in usable 
form, are found as the number of plates becomes 
infinite. Both of these steps are rather long and 
involved. A doubly infinite series must be set up 
and then summed. Series of Bessel functions of 
the zero and first order appear, which reduce 
eventually to the simple trigonometric terms of 
Eq. (12). 

It is of interest to note that the justifiable 
assumptions as to the magnitudes of the rotation 
and retardation found necessary in the previous 
method do not appear explicitly in this one 
though obviously they must be contained in it. 

From Eq. (12) we see that if the entering 
vibration coincides with an initial direction of 
principal stress, say the p direction, then the 
vibration along g at any point is given by 


R w 
w= (<1 sin ESx) cos («!- <x). (14) 
S VH 


If the ratio, R/2, of the rotation to the retarda- 
tion (due to the magnitude of stress) is small, 
even though the magnitude of the rotation itself 
be large, the intensity along g will be very small, 
the light will be approximately plane polarized 
at all points. The light vector may be said to 
follow along the direction of principal stress, i.e., 
practically no light will be scattered in this 
direction. The eye of an observer moving along 
a helix (the intersection of a line along the 
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Fic. 2. The light vector follows the axis of principal 
stress. O’g—entering direction of vibration of light coin- 
cides with a principal axis of stress. 


principal axis of stress with a cylinder having the 
points observed as its axis) will perceive a con- 
tinuous dark band. (See Fig. 2.) 

Equation (14) shows that even if the ratio of 
rotation to retardation is high, the intensity 
along g is zero wherever the vibration is plane 
polarized. There is a following around of the 
light vector, but between the points of plane 
polarization the intensity along g may be high. 
Fringes will therefore appear when there is 
rotation even though the entering vibration 
coincides with an initial direction of principal 
stress, while if there is no rotation, uniform 
intensity will prevail no matter what the direc- 
tion of observation. Furthermore, the spacing of 
the fringes will be the same as the spacing when 
the entering light has components along both 
p and q. 

The distance between fringes is decreased due 
to the rotation of the axes of principal stress. 
As can be seen from Eqs. (12) and (13), plane 
polarization or circular polarization will reoccur 
at intervals given by: 


(15) 


The retardation is thus seen to be increased by 
the factor S=(1+R?)! where R=2¢/A, or the 
distance between fringes decreased by 1/5. 
The difference between the secondary principal 
stresses in the general case is thus 

C7 

Cn" 
rf 


. (16) 
;S 
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which is to be compared with Eq. (1). As S can 
be quite large, neglect of it may, in many cases, 
produce appreciable error, though on the ‘“‘safe”’ 
side as the apparent stress will be higher than 
the true. If the total rotation is small, the 
formulas show that in practical cases it may be 
neglected as far as its effect on the distance 
between fringes is concerned. Usually one is 
interested in regions of high stress where the 
fringes will be closely spaced. In such cases, the 
total rotation probably will not be very large 
and the ratio of rotation to retardation probably 
small, though often appreciable in effect. 

While true extinction will occur only if the eye 
follows the directions of principal stress in a 
helix, it is nevertheless true that for a fixed direc- 
tion of observation, minima will occur at the 
same spacing as given in (15) and although the 
fringe pattern will not be as distinct it will yield 
the same results. The amount of rotation between 
two fringes may be obtained by finding the 
difference in the orientations for which each 
appears darkest. 


ANALYSIS OF GENERAL STATES OF STRESS 


Having the interpretation of fringe spacing in 
terms of stress magnitude and rotation of axes, 
it is possible to consider the procedures proposed 
for making the necessary optical measurements. 
Weller’s suggestion of universal rotation of the 
model while under load would require a very 
complicated loading device. Also, the stress at 
only one point at a time would be determined. 
The suggestion is therefore advanced that whole 
planes be investigated. 

The state of stress at a point is specified by 
six quantities, three normal stresses and three 
shearing stresses across the planes on which the 
normal stresses act. As the effect of an isotropic 
stress is only to change the initial index of re- 
fraction, it cannot be detected by a photoelastic 
pattern and therefore only the differences be- 
tween the normal stresses plus the three shearing 
stresses, or the difference between the principal 
stresses and the planes on which they act, can 
be found. Without a priori knowledge as to the 
state of stress, five independent measurements 
must be made to 
quantities. 


determine the unknown 
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Let a plane be chosen. Consider some line in 
the plane as a reference direction. If a narrow 
beam of circularly polarized light is passed in 
the plane so that the wave normal is first in the 
direction of the reference line, then at one angle 
to it and next at another, and photographs are 
taken of the light scattered at right angles to 
the wave normals, three independent measure- 
ments will have been made for each point in 
the plane. Call the reference line the x axis, the 
‘plane the xy plane. The fringe spacings obtained 
for three positions of the wave normal do not 
completely determine the state of stress. Two 
more independent measurements are required to 
get 0:—Gy, Gy—O2, Try, Tyz, Trz- Two positions of 
the wave normal in the xz plane other than along 
x will supply the additional information required. 
If, however, the orientations of the secondary 
principal stresses can be determined or are 
known, one need work only in the xy plane, the 
three fringe spacings and the orientations are 
sufficient. 

The general problem of finding the orientation 
of the secondary or primary principal stresses is 
rather complicated. If the rotation ratio is small 
at all points, then plane polarized light entering 
along a direction of principal stress will indicate 
the directions of principal stress all along the 
wave normal. Only when the entering light is 
vibrating in such a direction will a continuous 
dark band appear. If the ratio is large, no such 
reliable way is available. A continuous dark band 
will not exist, even the entering direction of the 
principal stresses cannot be determined. Dis- 
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tinctness of fringe patterns cannot be used as a 
criterion ; minimum distinctness does not neces- 
sarily occur when the observer is looking along a 
principal stress direction, and in general it would 
be difficult to determine such a minimum with 
any degree of exactness. 

Fortunately, in the majority of cases the ratio 
of rotation to retardation will be small and the 
directions of the principal stresses can probably 
be determined. 


ANALYSIS OF PLANE STRESS 


In plane stress, the magnitude and direction 
of a principal stress is known at every point, 
(o.-=0 acts perpendicular to the plane of the 
model). There are thus only three quantities to 
be determined, the principal stresses in the plane, 
71, o2, and their orientation. 

Consider rectangular axes x, y in the plate 
(Fig. 3). o,, the normal stress in the y direction, 
is determined from the fringe spacing when light 
is passed in the x direction, ¢: and o, when the 
wave normal is at two other known orientations 
with respect to the x axis. Taking these positions 
symmetrically about x, at some angle y to it, 
and calling the distance between fringes obtained 
along these directions d,,d;,d,, we have ¢,=C d,, 
o:=C/d:, «,=C/d,. The three normal -stresses 
determine the principal stresses and their orienta- 
tion in the plane. If y is taken as 60°, the for- 
mulas become rather simple. They are entirely 
analogous to those developed for the equiangular 
strain rosette.® 


wr 


of 





), (17) 


®R. D. Mindlin, ““The equiangular strain-rosette,”’ Civ. Eng. 8, 546 (1938). 
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Fic. 3. Analysis of plane stress. Fringe spacings obtained 


from three positions of the wave normal in the plane 
determine the principal stresses and their orientation. 


Similar formulas result for y= 45°. 

As o,+o,=0;1+02, two measurements at right 
angles are sufficient to determine the sum of the 
principal stresses. It will be more accurate for 
two-dimensional work to take an isochromatic 
picture in the usual manner and then find the 
sum of the principal stresses as above, using 
the additional result as a check on the accuracy 
of measurement. For the principal stresses them- 
selves, or for their sum and difference, care will 
have to be taken with the sign of the distance 
between fringes. If called positive when due to 
tension, it is negative when due to compression. 
The sign will often be apparent, if it is not, 
white light may be employed, the color sequence 
will show whether the stress is of the same or 
opposite sign as the stress at the previous point, 
and following back to a known point, as for 
example one on the boundary, the sign can be 
determined. Due to dispersion, it may be neces- 
sary to reduce the load for this determination. 
The isoclinics can also be used for determination 
of the sign. In the case of plane stress, the varia- 
tion of stress through the thickness is small and 
a wide beam of light can be employed. The wider 
beam will give a much greater intensity of 
scattered light and may therefore be desirable. 

It should be noted that a sheet of photoelastic 
material may appear free from initial bire- 
fringence when observed in a circular polariscope 
and yet have a fairly high equivalent stress 
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perpendicular to the plane, which may make it 
unsuitable for this proposed method. 

A two-dimensional state of stress also exists at 
an unloaded plane boundary. An isochromatic 
picture can be obtained to a certain extent but 
it would have little accuracy. The principal 
stresses and their orientation can be found by 
the method and formulas, (17), given for plane 
stress. 

A plane of symmetry or a principal plane in a 
three-dimensional case is very similar to the 
plane stress case. The quantities found from the 
three measurements at 60° will be the differences 
between the principal stresses and their orienta- 
tions, which is all that one can expect to obtain. 
Calling the normal stress on the plane az, just 
substitute o;—o- for o; and o2—<a; for a2 in (17). 


THE FREEZING METHOD 


The preceding development has all been 
applied specifically to the determination of the 
state of stress from the patterns produced by 
the light scattered from a plane. It applies also 
with but few modifications to the ‘‘freezing”’ 
method. In this method’* the model is heated, 
loaded, then cooled under load. Most of the 
birefringence present when the model is under 
load at the elevated temperature remains after 
cooling, and it is not disturbed by careful sawing. 
Thin slices are cut from the model and analyzed. 
When the slice is placed in a polariscope, an 
integrated effect through the thickness is ob- 
tained. If the stress in the neighborhood of a 
point is constant, the fringe order appearing 
there when the slice is in a circular polariscope is 


(dp—o,)t 
n= ; (18) 
( 

where, as before, o, and o, are the secondary 
principal stresses in the plane perpendicular to 
the wave normal (in general not the plane of the 
slice), tis the thickness of the slice in the direction 
of the wave normal, and C is the stress-optical 
coefficient. The difference between the secondary 
_ 7G, Oppel, ‘‘Polarisationsoptische Untersuchung raum- 
licher Spannungs- und Dehnungszustande,” Forschung auf 

dem Gebiete des Ingenieurwesens 7, 240 (1936). ; 
§M. Hetényi, ‘““The application of hardening resins in 


three-dimensional photoelastic studies,” J. App. Phys. 10, 
295 (1939). 
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principal stresses is therefore determined directly 
from the fringe order, their orientation can be 
obtained in the same way one obtains the iso- 
clinics for a model under plane stress. If the 
fringe orders and the orientations are obtained 
for three different directions of the wave normal, 
all information possible as to the state of stress 


is determined. The method is entirely analogous 
to that indicated in the section on ‘Analysis of 
General States of Stress’’ and applies with 
possibly even more accuracy because of the 
greater ease of determination of the directions of 
the secondary principal stresses. Little difficulty 
is caused by a small stress gradient provided the 
directions of ¢, and ¢, remain unchanged along 
the wave normal. 

Rotation of the slice about only one axis 
(can be taken in the slice) is required, and only 
three positions need be investigated in this 
rotation. Immersion in a liquid of approximately 
the same refractive index will eliminate re- 
flections. 

If the p and g axes do rotate, it should be 
possible to find the amount of this rotation in a 
plane polariscope if the stress is sufficiently high, 
and generally one is interested only in the highly 
stressed portions. As shown in the preceding 
work, plane polarized light entering along a 
principal stress direction will rotate with the 
principal stress axes and emerge practically plane 
polarized. Independent rotation of the polarizer 
and analyzer will find the orientations and the 
amount of rotation. The correction factor to 
the fringe order to find the stress is, as before, 
1/S. The difference between the secondary 
principal stresses is given by 

nC 1 


Tn—Oq>= ei (19) 


$ Ss 


The corrected values for the differences between 
the secondary principal stresses and the orienta- 
tions found as above should be used referred to 
the central plane of the slice. 

Hiltscher® proposed analyzing the slice cut out 
from the model in a microscope point by point 
using convergent polarized light as one does in 
crystallographic investigations. If rotation of the 
principal axes through the slice is at all appreci- 
able, the proposal does not appear practicable. 


CONCLUSION 


The investigation of the effect of rotation of 
the principal axes in planes perpendicular to the 
wave normal led to two important results: 

(1) Rotation increases the relative phase re- 
tardation by a factor S=(1+(2¢ A)?)! which at 
times may be quite large. If the ratio 2¢/A is 
not constant, the ratio of the increment in 2¢ to 
the increment in A between fringes can be used 
in regions of high stress. 

(2) At high stress, the light vector rotates with 
the principal stress directions, plane polarized 
light entering along a direction of principal stress 
will remain practically plane polarized in the 
direction of this stress through the entire model. 

The method of working in planes rather than 
point by point should be of value for general 
states of stress when a large portion of the model 
is to be analyzed. In the particular case of plane 
stress it yields a purely optical method of deter- 
mining the sum of the principal stresses, or the 
principal stresses themselves, over the entire 
model. 


® R. Hiltscher, ‘‘Polarisationsoptische Untersuchung des 
raumlichen Spannungszustandes im konvergenten Licht;” 
Forschung auf dem Gebiete des Ingenieurwesens 9, 91 
(1938). 
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Atomic Distribution in Aluminum-Silver Alloys During Aging 


CHARLES S. BARRETT AND ALFRED H. GEISLER* 
Carnegie Institule of Technology, Pittsburgh, Pennsylvania 
(Received June 17, 1940) 


Laue photographs of aluminum-rich aluminum-silver crystals during aging at 20° and 150°C 
show more or less clearly defined streaks along certain zonal ellipses, in addition to the streaks 
caused by thermal agitation. These have been analyzed by stereographic and reciprocal 
lattice projections and are found to result from two-dimensional gratings parallel to planes 
of the form {111}. The streaks with 20° aging are more diffuse than with 150° aging; no streaks 
were found with 200° aging. When the precipitate lattice is fully developed, it yields sharp 
spots along the streaks at positions predicted from earlier studies of the Widmanstatten 
structure; preceding this state there exist very thin plate-like nuclei on randomly spaced {111} 
planes, which govern the orientation and shape of the fully developed precipitate. These may 
consist of (a) clusters of silver atoms or (/) imperfect lattices caused by (111) layers of atoms 
shifting parallel to themselves as required for the transformation from the face-centered 
cubic to the hexagonal close-packed lattice but with the shifting occurring on random planes, 
thus destroying the lattice periodicity normal to these planes. 


INTRODUCTION 


NUMBER of instances have recently been 

discovered in which diffraction occurs with 
one of the three Laue conditions relaxed, the 
diffraction pattern resembling that of a two- 
dimensional, crossed grating. The phenomena 
accompanying precipitation in aluminum-rich 
aluminum-copper alloys are of this type. Guinier! 
and Preston? have shown that soon after an alloy 
is quenched from a temperature within the solid 
solution range a set of streaks begins to appear 
superimposed on the normal diffraction pattern. 
The streaks arise from two-dimensional gratings 
on planes of the form |100}—presumably from 
small plate-like clusters of copper atoms which 
are probably the nuclei of the @’ lattice. The 
clusters appear to be a few angstroms in thickness 
and several hundred angstroms in diameter and 
are dispersed at random through the matrix; as 
they lack periodicity in one dimension they give 
rise to diffraction maxima characteristic of two- 
dimensional gratings. As aging continues, the 
clusters become fully developed 6’ and ultimately 
the stable precipitate, 6. Some difficulties were 
encountered in the complete analysis of the 


* Aluminum Company of America Graduate Fellow. 

1 A. Guinier, Comptes rendus 204, 1115 (1937); 206, 1641 
(1938); Nature 142, 569 (1938). J. Calvet, P. Jacquet and 
A. Guinier, J. Inst. Metals 6, 177 (1939). 

2G. D. Preston, Nature 142, 570 (1938); Phil. Mag. 26, 
855 (1938); Proc. Roy. Soc. A167, 526 (1938); Proc. Phys. 
Soc. 52, 77 (1940). 
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photographs until it was realized that thermal 
agitation also causes streaks in Laue photographs 
(in the radial direction) which must be taken into 
account.* # 

Hendricks® has recently studied closely analo- 
gous streaks in x-ray photographs of various 
micas and has shown that they are produced by 
the lack of lattice periodicity normal to the 
cleavage plane of the mica crystals. Bradley, 
Bragg and Sykes® have observed similar diffrac- 
tion effects from submicroscopic lamellar arrange- 
ments of phases in iron-copper-nickel alloys with 
certain rates of cooling. Crossed-grating effects 
are also common in electron diffraction photo- 
graphs of very thin crystals and crystalline 
deposits. 

The following work is an analysis of Laue 
photographs of aluminum-rich aluminum-silver 
alloys during aging. It is shown that crossed- 
grating diffraction effects arise in the early stages 
of the precipitation process, somewhat as in the 
aluminum-copper alloys, but with important 
crystallographic differences. The crossed gratings 
are parallel to {111} planes of the solid solution 
in the present case, and are thus directly related 
to the habits of the precipitate crystals when they 








3G. D. Preston, Proc. Roy. Soc. 172, 116 (1939). 

*W. H. Zachariasen, Phys. Rev. 57, 597 (1940); S. 
Siegel and W. H. Zachariasen, Phys. Rev. 57, 795 (1940) 

5S. B. Hendricks, Phys. Rev. 57, 448 (1940). 

6A. J. Bradley, W. L. Bragg and C. Sykes, Researches 
into the Structure of Alloys, advance copy (Iron & Steel 
Inst. May, 1940). 
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F1G. 1. Projection of reciprocal lattice points and lines onto 
reference sphere centered at the origin. 


have grown to visible size, for Mehl and Barrett’ 
have shown that they are then plate-like and 
parallel to {111} planes. On the other hand, 
aluminum-copper, which contains crossed grat- 
ings on |100} planes in the early stages of aging, 
shows plate-shaped precipitate crystals on | 100! 
planes when they have grown to visible size.* 
Two different types of atomic arrangement that 
are capable of giving the effects are presented in 
the discussion. 


EXPERIMENTAL PROCEDURE 
Crystal preparation 


An aluminum-silver alloy of 20.2 percent silver 
prepared by melting in a graphite crucible under 
a potassium chloride flux was cast into an iron 
mold, rolled to a reduction of 50 percent in 
thickness, and homogenized 48 hours at 550°C in 
lampblack. Large grains were grown by the 
strain-anneal method, cut out with a jeweler’s 
saw, and thinned by careful rubbing on emery 
paper and by deep etching until the strained 
metal was removed. Crystals of chosen orientation 
were given a solution heat treatment of 18 hours 
at. 560°C in lampblack, quenched in iced water, 
and aged for various times at temperatures 
between room temperature and 320°C. 


X-ray photographs and their analysis 


Laue photographs were found convenient for 
the work and were made with the general radia- 
tion from a cobalt target tube having Lindemann 





7R. F. Mehl and C. S. Barrett, Trans. Am. Inst. Mining 
and Metallurgical Eng. 93, 78 (1931). 

§R. F. Mehl, C. S. Barrett and F. N. Rhines, Trans. 
Am. Inst. Mining and Metallurgical Eng. 99, 203 (1932). 
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windows and operating at 37 kv. The low in- 
tensity streaks required exposures roughly ten 


times as long as the normal Laue spots. 

For interpreting the Laue patterns, the recipro- 
cal lattice is useful’ and has been combined here 
with the stereographic projection which possesses 
some advantages over the usual gnomonic pro- 
jection for the purpose.'® Figure 1 is a sketch of 
the relations used. A reference sphere is centered 
at the origin of the reciprocal lattice, O, and the 
points of the reciprocal lattice are projected 
radially onto the sphere, as shown. These are 
projected, in turn, onto a stereographic projection 
by means of a projection point located on the 
surface of the sphere. The pole of a crystal plane 
on the reference sphere thus coincides with the 
radial projection of the reciprocal lattice point 
for that plane. 

If there is a lack of periodicity of scattering 
power in the lattice so that one of the three Laue 
conditions is relaxed and diffraction effects are 
those of a two-dimensional grating, the reciprocal 
lattice contains lines normal to the plane of the 
grating. In Fig. 1 this condition is indicated for a 
grating in the YZ plane. Heavy lines parallel to 
X, normal to the YZ plane, are drawn through 
the reciprocal lattice points (since it is assumed 
that ‘destruction of the third Laue condition 
leaves the others unaltered). The projection of 
these lines on the reference sphere is indicated in 
the figure, each point on the line being projected 
radially, as for the reciprocal lattice points. The 

















Fic. 2. Laue photograph of Al-Ag crystal aged 104 hours at 
150°C, beam along [110], with [001 ] vertical. 


*C. C. Murdock, Zeits. f. Krist. A99, 205 (1938). 
°C. S. Barrett, Trans. Am. Inst. Mining and Metal- 
lurgical Eng. 124, 29 (1937). 
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projected lines on the reference sphere (and on its 
stereographic projection) are zone circles passing 
through the point P, the pole of the crossed 
grating. 


RESULTS 


Figure 2 isa Laue photograph of an aluminum- 
silver crystal after quenching and aging 104 hours 
at 150°C. The x-ray beam was directed along the 
[110 ] direction. Each of the strong normal Laue 
spots is pierced by one or more streaks that lie 
along zonal ellipses. The geometry of these 
various streaks can be seen from the stereographic 
projection of the photograph, Fig. 3. (The method 
of plotting points along each streak is the same as 
if the points were strings of Laue spots.) It will 
be seen that all streaks follow great circles 
through the poles of planes of the type {111}, and 
correspond to the reciprocal lattice lines sketched 
in Fig. 4, which are parallel to directions of the 
form (111) in the reciprocal lattice; i.e., lines 
between (Akl) and (h4+1, R41, 1+1). Diffraction 
effects of this type would be produced by a set 
of crossed gratings arranged on the different 
}111} planes of the solid solution. 

Figure 5 is a Laue photograph with the x-ray 
beam along a cubic axis. Its projection in Fig. 6 
shows that the streaks follow zones through poles 





Fic. 3. Stereographic projection of one quadrant of Fig. 2 
using scheme of Fig. 1. Crossed-grating streaks (heavy 
lines) follow zone circles (dashed lines) through {111} poles. 
Radial asterism from thermal agitation is neglected. 
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Fic. 4. Sketch of reciprocal lattice for Figs. 2 and 3. 
Lines between (hk/) and (h41, k4+1, /41) are drawn as 
heavy lines where they correspond to the observed streaks 
of the Laue pattern. 


of {111} planes, as before. Photographs of other 
symmetrical and unsymmetrical orientations 
need not be reproduced here for they lead to the 
same conclusions. 

The early stages of aging at low temperatures 
usually produce a somewhat different appearance 
in the Laue photographs; for example, Fig. 7 
shows the appearance after fifty days at room 
temperature. A similar pattern was found in 
another crystal after one to four days aging at 
150°C. Four intense but diffuse maxima occur at 
the three, six, nine, and twelve o’clock positions 
near the center. Other streaks form a double 
cross near the center, the axes of the cross lying 
45° from the horizontal. These effects, however, 
can be accounted for on the same basis as the 
previous ones, if it is assumed that all streaks in 
the pattern are more diffuse than they are in 
Fig. 5. The greater width of the streaks through 
{131} spots causes them to overlap and blend 
with the thermal streaks at the three, six, nine, 
and twelve o'clock positions, and to diffuse into 
the background between these maxima and the 
{131} spots. The double cross can be ascribed to 
blending of several streaks which would corre- 
spond to great circles from (120) and (120) poles 
to the four {111} poles. 

On (110) Laue photographs the diffuse pattern 
contains radial rays extending out from the 
center and meeting the streaks through (121) and 
(133). In Fig. 4 these correspond to the dashed 
reciprocal lattice line from (000) toward (111). 
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Fic. 5. Laue photograph of Al-Ag crystal aged 104 hours at 
150°C, beam along [001], with [100] vertical. 





Fic. 6. Stereographic projection of one quadrant of Fig. 5. 


Now this line is (angent at (000) to the reflection 
spheres for all wave-lengths and therefore could 
diffract only in the forward direction were it not 
for the diffuseness of the reciprocal lattice lines, 
which permits contact with the spheres out for 
some distance from the center. The intermediate 
state giving diffuse streaks was observed with 
aging at 150°C and below, but not at higher 
temperatures. 

To show that Laue photographs are adequate 
to determine the plane of the two-dimensional 
grating, a Laue photograph of an aluminum- 
copper alloy is included, Fig. 8.* The alloy con- 
tained 4 percent copper and had been aged three 
days at room temperature after quenching; 


* This photograph was kindly given us by G. D. Preston, 
and is similar to, but more distinct than, some he has 
published. 
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radiation from a silver target was directed along 
[110 ]. The prominent streaks differ from those in 
Fig. 2, and their analysis by stereographic and 
reciprocal projections, Fig. 9 and Fig. 10, confirm 
Preston’s and Guinier’s conclusion that the two- 
dimensional grating is on cube planes in this 
system. In Fig. 10 the reflection sphere for 
Ag Ka is indicated; intersections of the sphere 
with the reciprocal lattice lines near (002), (024), 
(242) and (220) account for intensity maxima 
along the corresponding streaks in Fig. 8. 

In the aluminum-silver experiments, when 
aging is continued for about three hours at 200° 
or for shorter times at higher temperatures, new 
spots begin to appear along the Laue streaks, as 
in Fig. 11. Powder photographs of the alloy at 
this stage show that the normal precipitate has 
formed. (The structure of the precipitate is 
hexagonal close-packed, with a» = 2.879, c=4.573, 
and has a random distribution of silver and 
aluminum atoms.'') This is confirmed by the 
stereographic projection of Fig. 11, one quadrant 
of which is reproduced in Fig. 12. The new spots 
fall as predicted by the orientation relationships 
derived by Mehl and Barrett,’ namely, the (0001) 
plane of the precipitate parallel to |111} of the 
solid solution, with closest packed directions in 
these two conjugate planes parallel to each other. 

If the distribution of intensity along the streak 
through (012) and (133) in Figs. 2 and 3 is noted, 
it will be seen that the intensity drops to a 
minimum at 3 to § of the distance from one Laue 
spot to the next. The wave-length along this 
streak is almost constant, so the blackening is 
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Fic. 7. Laue photograph of Al-Ag crystal aged 50 days at 
room temperature, beam along [001], with [100] vertical. 


11 A. F. Westgren and A. J. Bradley, Phil. Mag. 6, 280 
(1928). ’ 
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here more nearly representative of scattering 
power than in the case of the other streaks. By 
the principles of diffraction from a grating con- 
taining a small number of lines, it follows that the 
position of the minimum corresponds to a thick- 
ness of the diffracting layers in the neighborhood 
of three to five unit distances. 

Powder photographs of polycrystalline samples 
quenched and aged for short times showed a few 
very weak extra lines not accounted for by either 
the solid solution or the precipitate. This indi- 
cates that a transition lattice! is playing a part 
in the mechanism, but the details are not yet 
clear. 


DISCUSSION AND CONCLUSIONS 


The evidence is conclusive for a type of 
crossed-grating diffraction from nuclei on the 
various {111} planes in the early stages of 
precipitation. The gratings appear to be between 
one and five unit cells thick with most heat 
treatments. Their lateral extent varies with the 
temperature of aging, and is of the order of 50 
unit cells at room temperature and many times 
greater than this at 150°C, as is indicated by the 
width of the streaks. 

Of what do the nuclei consist ? Two possibilities 
are presented here, and still others may be 
possible. The first is the analog of the theory 
presénted by Preston and by Guinier for alumi- 
num-copper alloys, and the second is a derivative 
of the mechanism postulated by Mehl and 
Barrett? for precipitation in the aluminum-silver 
system, and similar to one that has been sug- 
gested by Thomson and Cochrane™ to explain 
some anomalous effects in the electron diffraction 
from hexagonal close-packed metals that are 
transforming to their face-centered cubic modi- 
fication. 

The first theory is that silver atoms concen- 
trate in thin plate-like clusters that lie parallel to 
various {111} planes in different regions of a 
solid solution crystal and that are spaced ap- 
proximately at random. These clusters can be 
viewed as the first stage in the complex precipi- 


2 L. Guillet and L. Guillet, Jr., Comptes rendus 209, 79 
(1939). 


8G, P. Thomson and W. Cochrane, Theory and Practice 
of Electron Diffraction (Macmillan, 1939), p. 135. 
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Fic. 8. Laue photograph of Al-Cu, 4 percent Cu, aged 
three days at room temperature, beam along [110], with 
[001 } vertical. 





Fic. 9. Stereographic projection of one quadrant of Fig. 
8. Crossed-grating streaks follow zone circles through { 100} 
poles. 


tation process," identical with or followed by 
a stage in which a transition lattice, comparable 
to #’ in aluminum-copper alloys, is developed, and 
ultimately by the normal precipitate lattice. It 
has been argued" that these intermediate stages 
are the result of the complexities of atom move- 
ments by which the matrix lattice may transform 
to the equilibrium precipitate lattice under the 
crystallographic restraints imposed by the matrix 
lattice. As these alterations occur, the average 
size of the transforming regions is increasing 


“4 R, F. Mehl and L. K. Jetter, Preprint No. 30, October 
meeting (1939), American Society for Metals. 
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especially laterally. In any event, the basal plane 
(0001) of the precipitate can form from the 
(111) plane of the face-centered cubic matrix by 
merely shifting certain of the (111) layers of 
atoms parallel to themselves, and without altering 
the directions of the closest packed strings of 
atoms in these planes.’ In different regions of a 
crystal, this process occurs on different planes of 
the {111} type, and in each region a distinct 
orientation of the precipitated crystallite occurs, 
giving a total of four orientations for the precipi- 
tate. The precipitated particles retain their plate- 
like habit in growing toa visible size, and produce 
a Widmanstatten pattern of thin plates on {111} 
planes. This, in connection with the analogous 
behavior in aluminum-copper alloys* (in spite 
of the occurrence of a transition phase in the 
latter system) suggests that this is a general rule 
and lends new interest to microstructures of the 
Widmanstatten type. 

The second theory for the origin of the crossed- 
grating pattern is that the phase transformation 
starts at a great number of randomly distributed 
nuclei. In the early stages of the aging process, 
the alloy consists, on this view, of thin lamellae of 
the hexagonal close-packed phase distributed at 
random spacings through the face-centered lat- 
tice parallel to (111) planes of the latter. This 
could result from local regions in which a (111) 
net plane glides over the plane next to it as in 
crystallographic slip in such a way that the 
atoms in one plane are no longer directly above 
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Fic. 10. Sketch of reciprocal lattice for Figs. 8 and 9. 
Lines between (hkl) and (h+1, k,l), (h, k21,/) and 
(h, k, 141) are drawn as heavy lines where they correspond 
to the streaks of Fig. 8. 
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Fic. 11. Laue photograph of Al-Ag after precipitation (aged 
one hour at 320°C), beam along [100], and [001 ] vertical. 


atoms in the third underlying plane.’ The amount 
of the shift required for this is the distance be- 
tween neighboring hollows in the close-packed 
array of atoms on the plane. If shifts of this kind 
take place with proper regularity, they alter the 
sequence of planes of the face-centered lattice, 
ABCABC-: --,intothe sequence for the hexagonal 
close-packed, ABABAB:---. If, on the other 
hand, shifting is random, it leads to a destruction 
of periodicity, for atoms will be directly over 
each other at distances nd where d is the spacing 
of (111) planes and n=2, 3, 4, 5 etc., at different 
places in the crystal. This will produce diffraction 
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Fic. 12. Stereographic projection of Fig. 11. Crosses are 
calculated positions of poles of precipitate and open circles 
are observed positions; A refers to crystallites precipitated 
on (111) of matrix solid solution, B on (111), C on (111) 
and D on (111). Filled circles, triangles, and squares are 
poles of matrix solid solution. 
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effects such as we have observed, and analogous 
to those observed in the micas by Hendricks.°® 
With further aging, a consistent scheme will 
replace the haphazard one in the same way that 
long-distance order replaces local order in super- 
lattices, and sharp Laue spots from the precipitate 
lattice will appear along the streaks when some 
precipitate particles are large enough in all 
dimensions to give three-dimensional diffraction. 
When the precipitate has visible size, this process 


. would be expected to produce the microstructures 


Which of these two theories describes the pre- 
precipitation nuclei the more closely cannot be 
stated at present; nor can it be concluded that 
these are the only kinds of nuclei capable of ex- 
plaining the diffraction effects. But, nevertheless, 
the essential features of the process are clear. 
In the early stages of aging there are plate- 
like nuclei of extreme thinness on randomly 
spaced {111} planes, and these govern the 
orientation and shape of the Widmanstatten 
precipitate. The nuclei consist of—or develop 








ed and Widmanstatten orientations actually ob- into—a transition phase which later decomposes 
ul. served. into the stable precipitate. 
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1e IX. The Effect of the Pressure Variation of Viscosity on the Lubrication of Plane Sliders 
e, 

al M. Muskat AND H. H. EvINGER 

- Gulf Research & Development Company, Pittsburgh, Pennsylvania 

an (Received August 7, 1940) 

” The Reynolds theory is applied to the calculation of the lubrication properties of plane 

Ig sliders—thrust bearings—of infinite width provided with lubricants whose viscosities increase 

it exponentially with the pressure. The friction coefficient, minimum film thickness, and lubricant 

n flow were calculated both for fixed wedge angle and pivoted sliders. The effect of the viscosity 


variation with pressure is determined by the magnitude of the dimensionless product of the 
viscosity-pressure exponent and the bearing load per unit area. The analysis shows that for 
each choice of this product there will be a limiting position of the pivot line of the slider or 
of the equivalent Sommerfeld variable at which the film pressures and friction forces will 
become infinite, and beyond which it will be impossible to operate the slider. Moreover this 
product is shown to be limited by a maximum value equa! to 2, which means that the absolutes 
maximum load per unit area which can be carried by such slider systems is equal to twice the 
reciprocal of the viscosity-pressure exponent. Specific calculations on the friction properties 
of bearings operating with lubricants of different viscosity pressure exponents give curves of 
friction coefficient vs. load or Sommerfeld variable quite similar to those observed in practical 
tests. At high loads or low values of the Sommerfeld variable the friction coefficient curves 
split and follow the behavior generally interpreted in terms of oiliness and boundary lubrication 
phenomena. 


INTRODUCTION film viscosity as to change materially the 
N the last several years the question has been lubrication behavior of the bearings. It has been 
seriously raised as to whether or not the known for some time that the viscosities of 
variation of the viscosity of a lubricant with the ™ineral oils are increased upon the application 
fluid pressure may play a significant role in Of high fluid pressures, but that the rise in 
bearing lubrication in which the loads are viscosity due to this effect would be of im- 
abnormally high. In particular, it has been portance in practical lubrication problems has 
m" proposed that in the operation of slow speed and only recently been considered seriously. 


4 heavily loaded journal bearings, such as are used Among the first who have attempted to show 
1) in rolling mills, the high fluid pressures generated specifically that the rise in viscosity with 


re 


Ss 





in the lubricant film may so increase the normal 
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pressure can materially affect the performance of 
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lubrication systems are Bradford and Everett! 
and their collaborators. They showed that oils 
with different 
would give different operating temperatures and 
load carrying capacities when used in lubrication 


viscosity-pressure relationships 


systems under what were, otherwise, apparently 
identical conditions. Adding to their original 
interpretation, Givens*® found that the tempera- 
ture effects which they had observed with 
different oils could be directly related to the 
product of the viscosity-pressure coefficient and 
the load supported by the bearings. 

Another comprehensive study of this problem 
has been made by Needs.’ Needs measured 
directly the friction coefficient in a 120° centrally 
loaded fitted bearing as a function of the Som- 
merfeld variable‘ and load number of 
different lubricants, including mineral oil, castor 
oil, olive oil, and glycerine. In interpreting his 
experimental results, Needs assumed that his 
particular fitted journal bearing could be re- 
placed by an infinitely wide slider or thrust 
bearing. He then found that, whereas the theory 
for such sliders operating with lubricants of 
predicted 
decreasing monotonically with increasing load, 


for a 


fixed viscosity friction coefficients 
his data showed that as the applied load per 
unit area was increased the friction coefficients 
dropped to minimal values and then began to 
rise with continued increase of the load. The 
values of the load at which the minima were 
varied with the character of the 
lubricant. Needs then worked out the hydro- 
dynamic theoretical predictions for the variation 
of the friction coefficient with the Sommerfeld 


variable or bearing load of an infinitely wide 


observed 


slider for a lubricant whose viscosity increased 
exponentially with the pressure. The calculated 
friction curves for such systems were of exactly 
the same type as those he found experimentally. 
Although the quantitative agreement between 
the experiments and theory was not exceptionally 


‘L. J. Bradford and R. S. Wetmiller, Mach. Design, 
9, 36-39 (January, 1937); H. A. Everett, Soc. Auto. Eng. 
J. 41, 531, (1937); L. J. Bradford and C. G. Vandergrift, 
Inst. Mech. Eng. (General Discussion on Lubrication), 
23-29 (1937). 

2 J. W. Givens, Ind. Eng. Chem. 31, 1135 (1939). 

3S. J. Needs, Trans. Am. Soc. Mech. Eng. 60, 347 (1938). 

* The “Sommerfeld variable”’ is used here to denote the 
product of the lubricant viscosity and journal speed divided 
by the bearing load per unit projected area. 
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striking, the general features of the two were in 
complete agreement, and the orders of magnitude 
of the loads at which the minima in friction 
appeared were correctly given by the calculations. 
This concordance between the theoretical 
calculations and the experimental observations 
naturally appears to provide not only a con- 
firmation of the physical foundations underlying 
the hydrodynamic theory of lubrication, but 
also suggests that the increase of viscosity with 
pressure may, under certain conditions, provide 
at least a partial explanation for the commonly 
observed rapid rise of friction coefficient with 
decreased Sommerfeld variable, which is usually 
interpreted as resulting from the breakdown of 
viscous lubrication and the development of 
boundary lubrication phenomena. However, 
upon closer analysis of Needs’ theoretical con- 
siderations, it was found that the numerical 
calculations actually reported by Needs were 
restricted by severe limitations and assumptions. 
Specifically, all these calculations referred not 
only to a bearing 4” in length in the direction of 
sliding, but, moreover, it was explicitly assumed 
that the bearing is always so tilted that the film 
thickness at the leading edge is twice that at the 
trailing edge. As will be seen later in this paper 
such an assumption is not only artificial but is, 
moreover, impossible to sustain in the operation 
of any practical bearing if the viscosity does 
vary with the pressure. Because of this situation 
it seemed necessary that the problem be general- 
ized so as to make certain that even the qualita- 
tive agreement between theory and experiment 
so clearly exhibited in Needs’ work was not 
entirely spurious. On the other hand, Needs’ 
assumption that the sliders are of infinite width 
will be retained, since it has thus far not been 
possible to carry through all the necessary 
integrations for the finite width sliders. 
Although, as will be seen below, it is, in 
principle, possible by a simple transformation to 
take into account the variation of the viscosity 
with pressure as given by any type of analytical 
function or even by numerical or graphical 
representations of the variation, we shall actually 
assume, as did Needs, an exponential variation 
of the viscosity with pressure. The basis for this 
assumption lies in the fact that the direct 
empirical measurements of the viscosity-pressure 
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variation are generally best fitted by an expo- 
nential relation. Perhaps the closest correspond- 
ence between the empirical data and an expo- 
nential representation is that reported in the 
study of Hersey and Shore,* where the deviations 
of the semilogarithmic plots from linearity were 
generally rather small, and in some cases within 
the limits of the experimental errors. These data 
covered pressures to 4000 kg/cm? and tempera- 
tures to 140°C. 

In the present investigation it has, therefore, 
been explicitly assumed that the relation be- 
tween the viscosity and pressure is given by 
the equation: 

b= moe”, (1) 


where yy is the viscosity at atmospheric pressure, 


and ¢ is the vicosity-pressure exponent. The 
unit for the viscosity will be taken as the 
poise, and the pressure will be measured in 


kg/cm*. In these units, the published data give 
a range for ¢ of 5X10- to 3.5X10-* cm?/kg. 
In general, c decreases with increasing tempera- 
ture although the degree of variation with 
temperature is markedly different with different 
lubricants. The value of ¢ is in general greater 
for the mineral oils than for the fatty oils. 


ANALYTICAL THEORY FOR SLIDERS OF 
INFINITE WIDTH 


Using the notation of the seventh paper of 
this series® on sliders of finite width with fixed 
viscosity, we may write down at once the 
Reynolds equation for the infinite slider as 


d sh dp dh 
(— )=-0v : (2) 
dx\ uw dx dx 


5M. D. Hersey and H. Shore, Mech. Eng. 50, 221 (1928). 
Cf. also R. B. Dow, J. App. Phys. 8, 367 (1937), and M. D. 
Hersey, Theory of Lubrication ( (John Wiley, New York, 1938). 

6M. Muskat, F. Morgan, and M. W. Meres, J. App. 
Phys. 11, 208 (1940). 





VOLUME 11, NOVEMBER, 1940 





where (cf. Fig. 1) x is measured in the direction 
of motion, U is the velocity of the sliding plate, 
h is the film thickness, and uz is the lubricant 
viscosity, to be considered as a function of the 
film pressure p. 

Now, in principle, as previously indicated, 
any preassigned mode of variation of the vis- 
cosity with the pressure can be formally taken 
into account in the solution of Eq. (2) or the 
general Reynolds equation by simply introducing 
the new variable defined by 


Pp dp 
v= enn (3) 
n u(p) 


which reduces Eq. (2) to 


d dv dh 
—(#—)- —6U—. (4) 
dx dx dx 


This is identical in form with Eq. (2), if, in the 
latter, we suppose the viscosity to be fixed. 
Hence, if we simply solve the lubrication problem 
under the assumption of fixed viscosity the 
pressure distribution for the system in which 
the viscosity is variable will be given immediately 
by the inversion of Eq. (3), that is, by expressing 
p as a function of W according to Eq. (3). 

Here, however, we shall carry through the 
analysis explicitly only for a viscosity-pressure 
variation defined by Eq. (1), since, as we have 
already seen, the latter affords a reasonable 
approximation to the empirical data on the 
viscosity pressure relationship. The explicit 
transformation between WV and p as defined by 
Eq. (3) will, therefore, be given for our case by 


W=(1/"oc)(1—e-“”). (5) 

To solve Eq. (4) we now introduce the 

assumption that the film thickness / varies 
linearly with the distance x, that is, 

h=mx. (6) 


To fix the boundary conditions we note that 
since at the coordinates x; and xo, defining the 
limits of the lubricant film, the pressure (above 
atmospheric) must reduce to zero, the function 
WV must also vanish at these boundaries, in view 
of the definition of Eqs. (3) and (5). We thus 
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have at once for the function W the solution 


6’ (x; —x)(x—Xpo) 
y= (7) 
m*(x,+Xo) x? 


This solution, combined with Eq. 
directly the 


(5), gives 


the explicit form for pressure 
distribution in the film. 

To obtain the various derived properties of 
the system it is necessary only to carry through 
the appropriate integrations. First, however, it 


is convenient to introduce the following notation : 


0« pol l X; pol [ rf W 
B= ; p=-—; S= ; é= , (8) 
m?*(xo+x,) Xo Ww L 


where W is the total load per unit width sup- 
ported by the lubricant film and L is the slider 
length in the direction of motion, namely, x;— xo. 
The dimensionless function S once more appears 
in the analysis here as characterizing the com- 
posite dynamical reaction of the lubricant film, 
and is analogous to the similar dimensionless 
quantity appearing in the study of the journal 
bearings, and which has been termed the 
‘‘Sommerfeld variable.”’ The dimensionless quan- 
tity @ is a new variable introduced by the effect 
of the pressure variation of the viscosity, and 
its composition shows that for corresponding 
conditions of lubrication it is not merely the 
viscosity pressure exponent c which determines 
the behavior of the system but the product of c 
and the average film pressure W/L. 

By making use of Eqs. (5) and (7), and carry- 
ing out the integrations, it is found that the 
load supported by the film may be expressed as 


cW i ic¢ md 4p 
= f pdx =é= : 
L s To R+(p+1)" 
p+ 2/R pi 
x| log p—— tan '! | (9) 
p—1 p—1 VR 


where 6 has been eliminated by introducing the 
auxiliary variable R defined by 
R=4p B—(p—1)?. (10) 


The frictional forces may be obtained by 
integrating over the bounding surfaces the stress 
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components 


h ap 
o=pU/hF-—, 
2 Ox 


(11) 


where the upper sign refers to the moving plate 
and the lower one to the stationary slider. If 
out, and the 
fricticnal force be divided by the force normal to 
the plate, the coefficient of friction, f, for the 
moving slider is found to be given by 


f 1 Alf2pset+t (p+1)* 
=—+ ( ) tox p— | (12) 
nw 2 RL3#\p—-1 6 


In the case of the pivoted slider, it is necessary 


these integrations are carried 


to know the position of the center of pressure. 
This is found by imposing the condition that the 
summation of moments about the center of 
pressure vanishes. If Ax denotes the distance of 
the center of pressure from xo, it is given by 


Ax 1 “ 
= f (x — Xo) pdx 
L WL¥:z, 


1 2p 


~ (p—1)? &R+(p+1)2] 


*E(p?—1)(1+2@) —2p log p|— 
_ 


(13) 


Another interesting property of the slider is 
the lubricant flow between the sliders. This is 
given by the basic hydrodynamic expression for 
the liquid flow between parallel plates, and from 
which Eq. (2) was derived. For the pressure 


distribution appropriate to the case under 
consideration this is equal to 
Uh h® dp mUpL 
Q=—+ =——. (14) 


2 12udx = p*-!1 
The minimum film thickness, fo, is given 
simply by: 


ho=mxo=mL/(p—1). (15) 


An ‘auxiliary relation, useful in subsequent 
numerical calculations, may be obtained from 
Eqs. (8) and (10). By suitable substitutions 
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Fic. 2. The variation of f/m with S/m? for an infinitely 
wide slider operating with a lubricant whose viscosity 
increases exponentially with the pressure. f=friction 
coefficient ; m=slope of slider; S=poU/W; po=lubricant 
viscosity at atmospheric pressure; U=slider speed; 
W=load per unit width; ¢=cW/L; c=viscosity-pressure 
exponent; L=slider length, in the direction of motion. 
Dashed line is asymptotic limit of S/m?. 


making use of these relations it may be shown 


that 
| 2pf/pe+l 1 
simt=— (7) mneninorew, 
3@\ p—17 R+(p—1)? 


It should be noted that the quantity R was 
introduced merely as a convenience in writing 
some of the formulae. However, it turns out that 
its use also facilitates the numerical calculations. 

In order to obtain the numerical values shown 
graphically by the following curves, it was 
necessary to calculate sets of values of @, R, and 
p, which satisfy Eq. (9) and then to substitute 
these values into the formulae for the derived 
properties of the slider. Since @ was held constant 
for these curves, the most direct procedure was, 
of course, to substitute the desired value of @ 
into Eq. (9) and then to find values of R and p 
satisfying the equation. Newton’s method of 
approximation was employed for this purpose. 
For the curves presented in the next section in 
which the basic abscissa variable is S/m?, the 
calculations essentially consisted of an elimi- 
nation of p and R between Eqs. (9), (16), and one 
such as Eq. (12) or (14) for the various derived 
properties. For the section on the pivoted sliders 
the same type of elimination was carried out 


(16) 
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except that Eq. (13) rather than Eq. (16) was 
used to give the desired values of the abscissa 
variable. 


PLANE SLIDERS WITH FIXED WEDGE ANGLES 


As indicated by Eq. (16), the quantity S/m? 
enters the formulae determining the properties of 
a plane slider as a unit. However, it will be 
essentially equivalent to the directly observable 
quantity S itself if we consider the slider slope m 
to be fixed. From this point of view, the quantity 
f/m, as given by Eq. (12), may also be considered 
as representing directly the friction coefficient. In 
the next section, graphical representations of 
the properties of sliders will be given in which the 
sliders will be permitted to rotate about a pivot 
set along the center of pressure, and to assume 
various slopes depending upon conditions of 
operation. 

By the methods outlined above, calculations 
were made of the friction coefficient, given by 
Eq. (12), vs. S/m?. The results are plotted in 
Fig. 2 for various values of @. The curve ¢=0 
corresponds to the slider operating with a lubri- 
cant of fixed viscosity. It will be seen that with 
increasing values of @, that is, increasing viscosity 
pressure exponent, if the load be taken as fixed, 
the friction coefficient increases, although the 
general variation with S/m? for the larger values 
of the latter is very similar to that for the case 
é¢=0. This rise in the friction coefficient with 
increasing @ is clearly due to the rise in lubricant 
viscosity within the film resulting from the fluid 
pressures, which in turn require increased shear- 
ing forces in order to maintain the slider ve- 
locity U. 

At the low values of S/m? the curves for @>0 
show a radically different behavior than the 
curve for ¢=0. Thus, whereas the latter continues 
to decrease monotonically to the very limit, 
S/m?=0, the curves for @>0 fall to minimal 
values and then sharply rise again as S/m?* 
continues to decrease. In fact, the friction coefh- 
cients rise to infinity and become asymptotic to 
lines parallel to the friction axis at definite 
limiting values for S/m?. This peculiar behavior 
is the direct consequence of the general analysis 
of the previous section, and merely represents the 
fact that at these asymptotic limits the maximum 
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Fic. 3. The variation of the minimum film thickness /o 
with S/m? for an infinitely wide slider operating with a 
lubricant whose viscosity increases exponentially with the 
pressure. Notation as in Fig. 2. 


film pressures and associated film viscosities be- 
come infinitely high so that no further increase in 
load or decrease in S is permissible. This situation 
comes about analytically in the following manner. 
By combining Eqs. (5), (7), and (8) the pressure 
distribution in the film may be put in the form 
1 
e-? =—[ (1+8)y*— py(p+1)+Bp ], 
vy? 


(17) 


where y=x/xo. Hence, in order for the film 
pressure p to be real and have any physical 
significance, the bracket in Eq. (17) must be zero 
or positive. Recalling the definition of R from 
Eq. (10), it is readily verified that the bracket will 
always be positive for real values of y if R>O. 
Now at the limiting value of R, that is, when 
R=0, we have from Eqs. (9) and (16): 
4plogp SS  2p(p+1) 
t=—_——_; —=——__, (18) 
m* 3é@(p—1)? 
which, upon the elimination of p, give the relation 
between the minimal value of S/m? and é. Such 
an elimination shows, as can be inferred from 
Fig. 2, that the asymptotic limit for S/m?* 
monatonically increases with increasing 2, though 
it does not exceed 0.1! until @ becomes greater 
than 0.9. 
A very interesting corollary of the first of Eqs. 


(18) is that since p must always be real and 


positive—in fact greater than 1—é¢ can never 
exceed the value 2. The fact that the equation for 
é, as given by Eq. (18), refers to the limiting 
value for R, namely, R=0, gives the maximum 
possible value for @ under all conditions; for, by 
reference to Eq. (9) it will be seen at once that 
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when R>0O, the magnitude of @ will be even less, 
for any fixed value of p, than that given by Eq. 
(18). This limitation in the range of @ should be 
kept in mind in the consideration of the various 
curves, both of Fig. 2 and those of the later 
diagrams. 

The variation of the minimum film thickness 
with S/m? and for different values of @ is shown in 
Fig. 3. Here the rise in the values of ho/mLZ with 
increasing ¢ is clearly due to the fact that the 
increasing lubricant viscosity within the film 
makes it possible for the sliders to carry any given 
load or to operate with any fixed value of S/m? 
with a greater film thickness. In other words, the 
pressure increase in viscosity leads to an in- 
creased load carrying capacity for the lubricant 
film. 

Two curves for the variation of the center of 
pressure with S/m? are shown in Fig. 4. Curves 
for the other values of @ shown in the previous 
figures have been omitted because of the difficulty 
of keeping them separated on the scale used in 
the plotting. As a whole it appears from Fig. 4 
that the position of the center of pressure is 
relatively insensitive to the effect of the film 
pressure on the lubricant viscosity. To the extent 
that it has any effect it leads to a displacement 
of the center of pressure toward the trailing edge. 
The reason for this displacement is that the 
increase of the film pressure due to the viscosity 
variation is greatest in the region between the 
trailing edge and the center of the slider, thus 
leading to a resultant decrease in the moments 
of the pressure and with it in the distance of the 
center of pressure from the trailing edge. 
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Fic. 4. The variation of Ax/L with S/m? for an infinitely 
wide slider. Ax = distance of center of pressure from trailing 
edge. Rest of notation as in Fig. 2. 
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The variation of the lubricant fiow with S/m? 
is plotted in Fig. 5. The rise in lubricant flow with 
increasing ¢ is here caused by the increase in the 
film thickness already referred to. In fact, the 
effect of the larger film thickness is apparently 
more than enough to counterbalance the tend- 
ency for the lubricant flow to decrease because of 
the increased film viscosity. 


PIVOTED OR SELF-ALIGNING SLIDERS 


In most practical thrust bearing installations 
the sliders are pivoted so as to permit the slope of 
the lubricant film to adjust itself to whatever the 
conditions of operation may happen to be. In 
order for such an adjustment of the slope to take 
place, it is clear that the center of pressure in the 
lubricant film must fall underneath the pivot line. 
In other words, the pressure distribution within 
the film and its associated properties must adjust 
themselves so as to give a center of pressure 
coinciding with the pivot line. The position of the 
latter therefore essentially determines the detailed 
behavior of the lubricant film for any preassigned 
set of operating conditions. For this reason, in 
discussing the problem of the pivoted slider it is 
convenient to use the position of the pivot line 
as the principal defining parameter of the slider 
system. 

While it suffices to specify only the position of 
the pivot line to determine the lubrication 
properties of pivoted sliders operating with lubri- 
cants of fixed viscosity, it is necessary to assign 
also the value of @ in case the lubricant viscosity 
varies with the pressure. The various derived 
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Fic. 5. The variation of the lubricant flow with S/m?. 


Q=lubricant flow per unit width. Rest of notation as in 
Fig. 2. 
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Fic. 6. The variation of the slope, m, of an infinitely wide 
pivoted slider operating with a lubricant whose viscosity 
increases exponentially with the pressure. Ax/L = (distance 
of pivot line from trailing edge)/(length of slider in the 
direction of motion). Rest of notation as in Fig. 2. 


properties may, therefore, be plotted as a func- 
tion of the position of the pivot line for fixed 
values of 2, or conversely. 

The first derived property of interest for 
pivoted slider systems is that giving the slope 
assumed by the slider under various conditions of 
lubrication. It may be seen from Eqs. (16), (9), 
and (13) that this slope as well as the other 
derived properties of the lubrication system are 
essentially proportional to «/S, if any equivalent 
pairs of the other auxiliary variables are kept 
fixed. It is, therefore, sufficient to calculate or 
plot directly the coefficients of proportionality, or 
the ratios between the derived quantity and //S. 

The variation of m/./S with the position of the 
pivot line is shown in Fig. 6. A small variation of 
m//S with é@ will be evident. The general 
decrease of m/./S with increasing Ax/L is quite 
similar to that obtaining with sliders of fixed 
viscosity. In fact, the curve for €=0 has not been 
drawn on this figure because it would practically 
overlap that for ¢=0.1. The limitation in the 
range of Ax/L shown in Fig. 6 arises from the 
fact that the asymptotic limits of S/m? previously 
mentioned restrict the range of permissible Ax/L 
to continually increasing limits as ¢ is increased. 

The variation of the friction coefficient for the 
pivoted sliders may also be expressed in terms of 
the ratio f/./S, which is given essentially by the 
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Fic. 7. The variation of the friction coefficient, f, of an 
infinitely wide slider with the pivot position. Dashed line is 
asymptotic limit for pivot position. Rest of notation as in 
Fig. 2. 


product of the slope ratios of Fig. 6 and the ratio 
f/m of Fig. 2. The resulting curves plotted against 
Ax/L, with @ kept fixed, is given in Fig. 7. This 
set of curves shows that increasing the quantity é@ 
will result in increasing values of the friction 
coefficient. As may be inferred from the earlier 
discussion of the sliders with fixed wedge angle, 
this behavior is a consequence either of the 
increased film viscosity with increasing viscosity 
pressure exponent, if the average load is con- 
sidered as fixed, or the decreasing film thickness 
with increasing load if the viscosity pressure 
exponent be taken as fixed. As functions of Ax/L 
the f/.\/S curves all possess minima. However, 
with increasing ¢ the position of the pivot line for 
minimum friction coefficient, and hence more 
efficient operation, is shifted closer to the center 
of the slider. In the limiting conditions where 
¢@=2 the only permissible range of positions of the 
pivot line for stable operation apparently col- 
lapses and is shifted over to fall directly upon the 
center of the slider. 

The variation of the minimum film thickness 
with the position of the pivot line for fixed values 
of @ is shown in Fig. 8. With respect to the 
limitation of the ranges in Ax/L appearing in 
these curves, the same remarks apply as previ- 
ously made in connection with the function 
m//S. The general increase of the minimum film 
thickness for fixed S and Ax/L with increasing @ 
again is an expression of the fact that the in- 
creased lubricant viscosity due to the film pres- 
sures permits the film to sustain the applied load 
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with increasing film thickness, as would be 
expected from elementary considerations con- 
cerning slider systems operating with fixed 
viscosity lubricants. It will be noted from Fig. 8 
that the fy/Ly/S curves show maxima when 
plotted against Ax/L. Moreover, to the extent 
that the permissible range of Ax/L is large enough 
to allow for a full development of the maximum, 
the position of the latter shifts to decreasing 
values of Ax/L as @ is increased. Hence, in con- 
trast to the situation with respect to the friction 
coefficient, the effect of the pressure increase of 
the viscosity would here indicate that the pivot 
line be shifted toward the trailing edge in order 
to maintain maximum values of the film thick- 
ness and hence reduce the tendency for metallic 
contact and rupture of the lubricant film. 

The variation of the lubricant flow with Ax/L 
for fixed values of @ is shown in Fig. 9. Here, too, 
the curves as a whole show a trend which is in 
general parallel to that for the case ¢=0, in which 
the lubricant viscosity is independent of the 
pressure. To the extent that the permissible range 
of Ax/L permits a full development of maxima in 
the curves, the latter shift toward decreasing 
values of Ax/L with increasing ¢. The general 
increase of the lubricant flow with @ must here, 
too, be interpreted as being essentially due to the 
increasing film thickness with which the system 
can operate as a result of the increase of viscosity 
with pressure. 

In order to give a more practical representation 
of the significance of the various curves just 
presented, some calculations were carried through 
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Fic. 8. The variation of the minimum film thickness, ho, 


of an infinitely wide pivoted slider with the pivot position. 
Rest of notation as in Fig. 2. 
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on the properties of several slider systems 
operating under specific conditions. In all cases 
the atmospheric lubricant viscosity was taken as 
1 poise, the slider velocity as 10 cm/sec., and the 
pivot positions, Ax/L, as 0.4. Calculations were 
then made for the variation of the friction 
coefficient and film thickness for three lubricants 
having viscosity pressure exponents of 0, 10-4, 
and 2-10 in.?/Ib. Since the atmospheric vis- 
cosity and slider speed were fixed, the only 
physical variation which would give a variation 
in S was that in the average load. In Fig. 10 the 
resulting friction coefficient as a function of the 
average load is plotted in terms of the composite 
variable S. In order to avoid having to specify 
the absolute value of the slider length in the 
direction of motion the products f,\/Z and «/SL 
were used as units. From Fig. 10 we see at once 
that, whereas for the lubricant whose viscosity 
is independent of the pressure the friction coeffi- 
cient curve is a straight line passing through the 
origin, those for the lubricants with variable 
viscosity fall to a minimum and then rise steeply 
to infinity as S approaches its asymptotic per- 
missible limit. This result is clearly very similar 
to, and suggestive of, the almost universal 
empirical observations with all types of lubri- 
cation systems, namely that as the characteristic 
variable corresponding to S continues to decrease 
the friction coefficients ultimately reach mini- 
mum values and then sharply rise to values which 
may be enormously greater than those at the 
minimum. It is generally considered and accepted 
that the passage through the minimum and 
subsequent sharp rise reflects a development of at 
least partial metallic contact and the incidence of 
boundary lubrication conditions. Here—as has 
also been pointed out by Needs—we find an 
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; Fic. 9, The variation of the lubricant flow, Q, of an 
infinitely wide pivoted slider with the pivot position. Rest 
of notation as in Fig. 2. 
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Fic. 10. The variation of the friction coefficient with S 
of pivoted sliders operating with three kinds of lubricant. 
Atmospheric viscosity=1 poise for all lubricants; slider 
speed =10 cm/sec. and pivot line distance =0.4 of slider 
length in all cases; c=viscosity pressure exponent ; dashed 
lines indicate asymptotes limiting permissible values of S. 
Rest of notation as in Fig. 2. 


exactly similar type of behavior under full 
viscous hydrodynamic conditions without having 
to invoke the complicating factors obtaining 
under boundary lubrication conditions. It will be 
seen, too, that the exact position of the minimum 
and the numerical values of the friction coeffi- 
cient in the sharply rising portion beyond the 
minimum markedly depend upon the absolute 
values of the viscosity pressure exponent c, and 
in this way give a purely hydrodynamic analog of 
the so-called property of ‘‘oiliness.’’ In particular, 
Fig. 10 shows that the lubricants with smaller 
viscosity pressure exponents will be more “‘oily’’ 
than those with high exponents, if one were to 
interpret the curves in terms of the oiliness 
phraseology. Of course, it is not proposed that in 
general the hydrodynamic effect of the increase 
in viscosity with pressure gives a real explanation 
of what are commonly termed boundary lubri- 
cation conditions or of the oiliness properties of 
the lubricants. However, in such systems where 
the loads may be exceptionally high, it may well 
be that the passage of the friction coefficient 
curve through a minimum and subsequent rise 
may be due, at least in part, to the purely 
hydrodynamic effect discussed here. 
Corresponding to the higher friction coefficients 
indicated in Fig. 10 for the cases c>0, the mini- 
mum film thicknesses will also be greater. In fact 
whereas the ho vs. \/S curve will pass linearly 
through the origin for c=0, the film thickness for 
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Fic. 11. The pressure distributions for the pivoted slider 
systems described by Fig. 10. For curves on left, average 
load = 9400 Ib. /in.?; for curves on right, average load = 5000 
Ib./in.*; for solid curves c=2X10~™* in.?/Ib.; for dotted 
curves ¢c=0 (const. viscosity); (x—xe)/L = (distance from 
trailing edge) /(slider length, in direction of motion). 


c>0 will always remain nonvanishing up to the 
asymptotic limits of S. Hence, while the sharp 
rise in the friction coefficient on passing the 
minimum would lead to high energy losses in the 
continued operation of the system, the fact that 
the film thickness remains relatively large and 
nonvanishing should give a compensating effect 
in avoiding the scoring and scuffing of the sliding 
surfaces. 

The pressure distributions in systems corre- 
sponding to those represented by Fig. 10 are 
shown in Fig. 11 for c=0 and c=2-10~ in.?/Ib. 
For the curves on the left the average load is 9400 
lb./sq. in., which is the maximum permissibie for 
the lubricant of variable viscosity ; the curves on 
the right correspond to an average load of 5000 
lb./sq. in. It will be seen that the effect of the 
pressure increase of viscosity is to concentrate 
the normal film pressure rise to the region in the 
vicinity of the maximum obtaining in the film of 
constant viscosity. In fact, as indicated by the 
solid curve on the left, where the load is the 
maximum permissible, the film pressure theo- 
retically rises to infinitely high values. At the low 
loads, of course, the whole effect is not so marked 
although it is still quite appreciable. 

While Needs’ theoretical work was pretty 
much limited to a presentation of curves of the 
type shown here in Figs. 10 and 11 for specific 
numerical cases, the corresponding results derived 


here are as a whole quite similar to those given by 
Needs. His explanation of his experimental obser- 
vations, which show the same types of minima 
exhibited in Fig. 10, as being due to the effect of 
the pressure on the lubricant viscosity can 
therefore be accepted as sound. However, the 
present work does show that his specific assump- 
tion that the parameter p always has the value 2 
is not justified. In particular, for the cases 
illustrated in Fig. 10 the value of p actually varies 
from a minimum of 1.0 to a maximum of 1.85. 
On the other hand, as indicated at the beginning 
of this paper, even the work presented here is 
based on the assumption that the sliders are of 
infinite width, whereas in practice they must, of 
course, always be of finite width. 

It may be noted finally that the correlation 
observed by Givens between the temperature 
rise in the lubrication systems of Bradford and 
Everett and the product of the viscosity pressure 
coefficient and the load is quite consistent with 
the results of the present investigation. Although 
Needs did not point out the significance of this 
product, which has here been denoted by é@, the 
analysis of this paper shows that this quantity 
appears naturally as one of the defining parame- 
ters of the lubrication system. While this paper 
has not dealt with the temperature effects associ- 
ated with the lubrication problem it is clear that 
these effects, too, will be ultimately determined 
by the composite quantity é@ rather than the load 
or viscosity pressure exponent separately. On the 
other hand, it should be noted that these same 
temperature effects within the lubricant film will 
tend to lower the effective values of the exponent 
c and thus extend the observed limiting load 
carrying capacities beyond those predicted on the 
assumption that c remains fixed. Moreover, the 
general lowering of the viscosity due to the 
internal film temperature elevations also will 
undoubtedly modify the pressure distributions 
and resultant integrated properties of the prac- 
tical lubrication systems. 
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